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<I>-variational stability for discontinuous system 

Xuefeng Liang! and Wansheng He^ 

School of Mathematics and Statistics Tian shui 
Normal University, Tian shui 741000, Gan su, China 
E-mail: liangxf0113@126.com 

Abstract Bounded ® -variation functions are development and generalization of bounded 
variation functions in the usual sense. The concept of Henstock-Kurzweil integral is an ef- 
fective tool in dealing with highly infinite oscillation functions. In this paper, the concept of 
4>- variational stability is defined by using ® -function theory, the $-variational stability of the 
bounded <S>-variation solution to discontinuous system is discussed, and the Ljapunov type 
theorems for 4>-variational stability and asymptotically <I>-variational stability of the bounded 
4>-variation solution are established. These results are essential generalization of variation 
stability of bounded variation solution of this system, and the certain foundation is laid in the 
research of highly infinite oscillation functions. 

Keywords Discontinuous system, Henstock-Kurzweil integral, bounded 4>-variation solution, 
4>- variational stability, Ljapunov function. 

2000 Mathematics Subject Classification: 34A20, 34G10. 
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§1. Introduction 

Consider the general discontinuous system 

x' = (*) 

where x = (x\, X 2 , x n ) T , x' = ^f,/ : G — » R n , G is an open region in R n+1 . (*) is called 
discontinuous system if / is a function with some discontinuity. In [1] and [2], the existence, 
uniqueness and stability for the solution of Caratheodory system and Filippov system were 
obtained by using the Lebesgue integral and the solutions obtained are absolute continuous 
functions. The Henstock-Kurzweil integral was introduced by Henstock and Kurzweil indepen- 
dently during 1957-1958, the existence, uniqueness and variational stability for the bounded 
variation solutions of a class of discontinuous system were discussed by using the Henstock- 
Kurzweil integral in [3,4]. 

The functions of bounded ^-variation were introduced by Musiclak and Orlicz I 5,6 !, which 
are generalization and development of functions of bounded variation in usual sense. The 

1 This work is supported by the Young Teachers’ Scientific Research Projects of Tianshui normal university 
(TSY 201202). 
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existence theorem of bounded ^-variation solutions for a class of discontinuous system is es- 
tablished in [ 7 ]. In this paper, the concept of <f>- variational stability of bounded (^-variation 
solutions for discontinuous system is given and the stability is discussed through making use 
of Ljapunov function. The Ljapunov type theorems for $-variational stability and asymptoti- 
cally ^-variational stability of the bounded ^-variation solutions are established. These results 
are essential generalization of variation stability of bounded variation solution to discontinuous 
system in [ 4 ]. 

§2. Preliminaries and definitions 

Definition 2.lJ 3,4,7 >8] A function x : [a, 6] — » R n is called Henstock-Kurzweil integrable 
over [a, b], if there is a A G R n such that given e > 0 , there is a positive function S : [a, b] — » 
( 0 , Too) such that for any S(t)— fine partition D : a = t$ < t\ <■■■< tk = b and {6, £2, ■ ■ ■ , 6} 
satisfying 6 G [tj-i,tj] C [£,• - < 5 ( 6)>6 T < 5 ( 6 )]> we have 

k 

iiE^x^- 1 )-^ < e - 

3 = 1 

A is called the Henstock-Kurzweil integral of x(t) over [a, b] and will be denoted by f ' x(t)dt. 

Definition 2.2.I 3 ’ 4 ’ 7,8 ! A function x : [a, b] — > R n is called H-K-Stieltjes integral over [a, 6], 
function y : [a, b\ — > R n , if there is a A G R n such that given e > 0 , there is a positive function 
S : [a, b] — » ( 0 , Too) such that for any 6 (t)— fine partition D : a = <0 < f 1 < • • • < tk = b and 
{6,6, • • ' ,6} satisfying 6 G [tj-iAj] c [6 ~ + 66)]> we have 

k 

II ’ 52 x{Zj){y{tj) - y(tj- 1)) - A|| < e. 

'M 1 

A is called the H-K-Stieltjes integral of x(t) over [a, 6] and will be denoted by x(t)dy(t). 

Let $(it) denote a continuous and increasing function defined for u > 0 with $( 0 ) = 0 , 
4 >(u) > 0 for u > 0, and satisfying the following conditions: 

(Ai) There exist u 0 > 0 and a > 0 such that 4 >( 2 u) < o 4 >(u) for 0 < u < u 0 ] 

(A 2 ) $(u) is a convex function. 

Let [a, b] C R,— 00 < a < b < Too. We consider the function x : [a, b] R n , x{t) is of 
bounded ^-variation over [a, b] if for any partition n : a = to < t± < • • • < t m = b, we have 
V$(x; [a, b ]) = sup w Y^Li 'KII^’fA) ~ x(i*-i)||) < Too, V<s>(x; [a, b ]) is called ^-variation of x(t) 
over [a, b\. We always assume $(w) satisfying (Aj) and (A 2 ). 

Given c > 0 , we denote B c = {x £ R n ; ||x|| < c}. Let / C (a, b) C R be an interval with 
—00 < a < b < Too, and set G = B c x I. 

Definition 2.3. 1 3,4 , 7 ! A Caratheodory function / : G — > R n belongs to the class V$(G, h, w), 
if 

( 1 ) There exists a positive function S : [a, b] — ► ( 0 , Too), satisfying r G [u, v] C [r — S(t), rT 
<5 (r)] C I and x G Hfor every [it, u], we have 



|| f{x,r)(v-u) ||< ®(\h(v) - h(u)\). 
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(2) Satisfying r € [it, i>] C [r — 5(r), r + <5(r)] C I and x,y £ B for every [it, i>], we have 

II f(x,r) - f(y,r) || (v - u) < w(|| x-y ||)$(|/i(v) - h(u)\), 

where h : [a, b] — > R is an increasing function and continuous from the left on the interval [a, b ], 
u> : [0, +oo) — > R is a continuous and increasing function with w(0) = 0. 

(3) The function f(ip(t),t) is H — K integrable on [a,0\ for every step function if>(t) on 
[a,/3] C [a, b\. 

Lemma 2.1. [ 3 ’ 4 > 7 1 Assume that / : G — > R n satisfies the condition (1) of Definition 
2.3, if x : [a,0\ — > R n ,[a,/3\ C (a, b) is such that (x(t),t) £ G for every t £ [a,0\ and 
if the Kurzweil integral J 13 f(x(t),t) exists, then for every pair S\,S 2 € [a,/3] the inequality 
II la f(x{t),t)dt ||< V$(/i;[si,s 2 ]) holds. 

Corollary 2 .lJ 3 ’ 4 ’'] Assume that / : G — » R n satisfies the condition (1) of Definition 2.3, 
if x : [a,/3] — > R n , [a,/3] C (a, b) is a solution of (*), then x is of bounded ^-variation and 
V$(a;; [a,/?]) < $(V$(h; [«,/?])) < +oo. Moreover every point in [a,/3] at which the function h 
is continuous from the left is a left continuity point of the solution x : \a,0\ —> R n . 

Lemma 2 .2.I 4,8 ! Assume that — oo < a<b < +oo and that f,g : [a, b\ — > R are functions 
which are continuous from the left in [a,b\. If for every a £ [a, b] there exists 5(a) > 0 such 
that for every 77 £ ( 0 ,d(cr)) the inequality f(a + 17 ) — f(a) < g(a + 17 ) — g(a) holds, then 
f(s) - f(a) < g(s) - g(a) for all s £ [a, b\. 

If f(0,t) = 0 for every t, £ [0, + 00 ) then the function x given by x(t) = 0 for t > 0 is a 
solution of (*) on the whole half-axis [ 0 , + 00 ). 

Definition 2.4. The solution x = 0 of (*) is called <f>-variationally stable if for every 
e > 0 there exists 5 = d(e) > 0 such that if y : [to,ii] —> B c , 0 < to < ti < +00 is a function 
of bounded 4>-variation on [t 0 ,ti], continuous from the left on (f 0 ,ii] with $(||i/(to)||) < 5 and 
V<s>((y(s) - f t S 0 f(y(t),t)dt); [t 0 , ti]) < 5 , then we have $(||y(t)||) < e. 

Definition 2.5. The solution x = 0 of (*) is called <f>-variationally attracting if there 
exists So > 0 and for every £ > 0 there is a T = T(e) > 0 and 7 = j(e) > 0 such that if 
y : [to, ti\ —> B Cl 0 < <0 < < +00 is a function of bounded ^-variation on [to, ti], continuous 

from the left on (t 0 , ti] with $(||y(to)||) < <5o and V<s>((y(s) - f(y(t),t)dt ); [t 0 , ti]) < 7 , then 

<f > (||i/(t)||) < £ for all t £ [to, ti] f|[t 0 +T(£),+oo) and to > 0 . 

Definition 2.6. The solution x = 0 of (*) is called <f>-variationally asymptotically stable 
if it is <f>-variationally stable and <f>-variationally attracting. 



§3. Main results 

Theorem 3.1. Assume that V : [0, + 00 ) x R n — > R is such that for every x £ R n the 
function V(-,x) : [0, + 00 ) — > R is continuous from the left in (0, + 00 ) and 
(1) The inequality 

\V(t,x)-V(t,y)\<L$(\\x-y\\) (1) 



holds for x, y £ R n , t £ [0, + 00 ) with a constant L > 0; 
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( 2 ) There exists a real function H : R n — » R such that for every solution x : (a, / 3 ) — > R n 
of the discontinuous system (*) on (a,/ 3 ) C [ 0 , +oo), we have 



lim sup 
v^o+ 



V(t + V, x(t + 77)) - V (t, x(t)) 
V 



< H(x(t)) 



( 2 ) 



for t € ( a , ( 3 ). 

If y : [to, ti] — > R n , 0 < to < ti < +oo is continuous from the left on (to, ti] and of bounded 
^-variation on [to,ti], then the inequality 

V{ti,y(ti)) < V(t 0 ,y(t 0 )) + NV^((y(s) - [ f(y(t),t)dt);[t 0 ,t 1 ]) + R(t 1 -t 0 ) ( 3 ) 

Jt 0 



holds, where R = sup te [ to tl ] H(x(t)),N = Lx(a'), (x( a ') was defined by Theorem 1.2 in [ 5 ]. 

Proof. Let y : [t 0 , ti] — > R n , 0 < to < ti < +oo be a function of bounded ^-variation on 
[to, ti] and continuous from the left on (to, ti] and let c £ [to, ti] be an arbitrary point. It is 
clear that the function V(t,y(t)) : [to, ti] x R n — ► R is continuous from the left on (to, ti]- 

Assume that x : [c, c + 771 (o')] — > R n is a solution of bounded <!>- variation of (*) on the 
interval [c, c + 771(17)], 771 > 0 with the initial condition x(a) = y(c 7). The existence of such a 
solution is guaranteed by Theorem 3.3 in [ 7 ]. By ( 1 ), we have 



V(a + r),y{a + 77)) - P(c + 77, x(a + 77)) < L$(\\y(a + 77) - x(cr + 7?)||) 

r cr+ri 



= L${\\y{(T + 77) - y{a) - / f(x(t),t)dt\\) 

J G 



for every 77 £ [0, 771(17)]. 
By ( 2 ), we obtain 



V(er + 77, y(cr + 77)) — V (c, x(a)) 



= F(c + 77, y{cr + 77)) — V(a + 77, x(a + 77)) + V(a + 77, x(a + 77)) - V ( c , x(a)) 

< L<$>(\\y(a + t?) - y(a) - f° +r> f(x(t),t)dt\\) + T)H{x{a)) 

< L${\\y{a + t?) - y(a) - f° +r> f{x(t),t)dt\\) + r)R + t?£, 

where e > 0 is arbitrary and 77 € (0, 772(a)) with 772(c) < 771(c), 772(c) > 0 is sufficiently small. 

Denote s 

'&{s) = y(s)~ f f(y(t),t)dt,s £ [t 0 ,ti], 

Jt 0 

the function T : [t 0 ,ti] — > R n is of bounded ^-variation and continuous from the left on (t 0 , ti]- 
By Theorem 1 . 2 , Theorem 1.11 and Theorem 1.17 in [ 5 ], we have 

V(a + 77, y{a + rf))-V (c, x{a)) 

< L$(\\y{a + 77) - 1 /(c) - f° +ri f(y(t), t)dt\\ + || J^ + ' n [f{y{t),t) - f(x(t),t)]dt\\) 

+yR + 7 ?e < Lx{a'){§{\\y{<j + 77) - y{a) - f° +v f{y{t),t)dt \\ ) 

+ $ (ll XT +, '[/(t/W ,0 - /(x(t),t)]dt||)} + r]R + r]e 
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< Lx(a')V<s,{{y{s) - £ f(y(t), t)dt)\ [a, a + y}) + yR + ye 

+Lx{a ')§{ || S° +r, [f{y{t),t) - f(x(t),t)\dt\\) 

< VVb(^ ; [a,a + y]) +yR + ye + N<f>{\\ J° +V [f(y(t), t) - f(x{t),t)]dt\\) 

< 7V(V$(^; [,a + y]) - Vb(^; [,<j])) + yR + ye 

+Nm\C +ri [f(y(t),t)-f(x(t),t)}dt\\) (4) 

for every y £ ( 0 , 772 (cr)). 

Since / € V$(G, h,ui) and by Definition 2.1, there exists a positive function d(r) : [cr, a + 
77 ] — > (0, +00) for e > 0 such that any S(t)— fine partition D = {(tj; \otj- 1 , otj]),j = 1, 2 , • • • , to} 
on [<7, a + 77 ] , we have 

II Ja +V if{y(t),t) ~ f{x(t),t)]dt\\ 

< II J° +V [{y{t),t) - f(x(t),t)]dt - - f( x ( T j), a j) 

~f{y( x j)^j-i) + f( x ( x j),otj-i)]\\ 



+11 Ej=i[/(y(D)- a i) - - f(y( T i)’ a j-i ) + f( x (Tj)i a j-i) 



< §+£7=i“(ll ll)*(l I)- (5) 

By Theorem 1.3 and Theorem 1.17 in [5], we have 

£7=1 “(II ll)*(l I) 

= £7=1 mi i/(d) - x ( t j) «j]) 

< £7=1 “(II y( T i) - ^(D') ll)[Vb(*.; [o-.aj-]) - V${h; [a, ac^i])]. 

Denote U{t) = V$(/i; [cr, t] ) for cr < t < cr+y. By definitions of $(tt) and h, U(t ) is a function 
of non- navigate and nondecreasing on [cr, a + 77] and continuous from the left on (cr, cr + 77 ] , then 
by (5) we have 

£7=1 mil y(rj) xirj) ||)$(| h(aj) - h(a 3 . 1) |) 

< £7=1 “(II v(tj) x ( x j) I m(aj) - Uiaj.!)) 

<11 £7=1 “(II vfr) - *fo) ll)(C/K) - U(aj- 1 )) - w(|l 7/(t) - *(i) ||)dtf(i) || 

+ r + "“(ll y(t)- x (t) W)dU(t) < |+/7 +,, w(|| y{t)-x{t) II )dU(t). ( 6 ) 

By (5), ( 6 ), Theorem 1.16 in [ 8 ] and arbitrary e > 0, we have 

II J° +ri [f(y(t),t) - f(x(t),t)}dt\\ < /7 +,? w( II y(t)-x(t) W)dU(t) 

= lim a ^ 0 +[/7 + ““(ll y{t) - x (t) || )dU(t) + X77 w(|| y{t) - x(t) \\)dU(t)\ 

= “(II y(v)-x(a) ||)(G(ct+) - U{pj) + lim a ^o + /77 ? “(|| y(t) - x(t) \\)dU(t) 
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= lim Q ,_ J . 0 + w(|| y(t) - x(t) || )dU(t) 

< sup pe[(7)<T+77 ] w(|| y{p)-x(p) ||)lim a ^o+(f/(o- + ry) - U(a + a)) 

= suPpel^-H,] w(|| y{p)-x{p) \\)(U(a + r])-U(<j+)). ( 7 ) 

Let p G [a, a + 772(17)], we have y(p) — x(p) = y(p) — y(a) — ff f(x(t),t)dt , and therefore 
lim (y(p) - x(p)) = y(a+) - y(a) - lim ( f(x{a),p ) - f(x(a), a)) 

p— ><T+ p— ><T + 

= y(?+) - y(o-) - (f(x(a),a+) - f(x(a),c r)) = ^(cr+) - $ (a), 

and also 

lim || y{p) - x{p ) || = ||tE , (cr+ ) - 'F(o")||. (8) 

p->< T + 

For every e > 0 denote 

^ = N(U(h) - U(t 0 ) + 1 ) > ° (9) 

and assume that r = r(/ 3 ) > 0 is such that w(r) < f 3 . Further, we choice a G ( 0 , |). 

Since (8) holds, there is an 773(0") € (0,772(0")) such that 

I \y(p) - a;(p)ll < ll^(o'+) - ^WII + a, 

and also 

w(H y(p) - x(p)ll) < w(||W(cr+) - 'F(0")|| + a) (10) 



for p G (cr, a + 773(0-)). 

Denote P(/ 3 ) = {a G [to,ti]; ||\ 1 /(ct+) — \I/ (cr) || > |}, since is of bounded ^-variation on 
[fo,ti], the set P( 0 ) is finite and we denote by n(/ 3 ) the number of elements of P(/ 3 ). 

If o" G [to, ti] \ P{ft) and p G (a,a + 773(0")) then by ( 10 ) we have 

w(|| y(p) ~ *(P)II) < w(||^(ct+) - 'F(o")|| + a) < uiC- + a) < u>C- + *-) = w(r) < /?, 

and by ( 7 ) also 



rv+V 

11 / [f(y(t),t) - f(x(t),t)\dt ||< / 3 (U(a + r]) - U(a+)) ( 11 ) 

J a 

whenever 77 G (0,773(0")). 

If a G [fo,ti]n-P(/ 3 ) then there exists 774(0") G (0,773(0")) such that for 77 G (0,774(0")) we 

have 



U(cr + 77) — U(a+) =| U(a + 77) — U(cr+) \< 



P 



(n(P) + l)w(||^(o"+) - ^(ct)|| + a) 
and (a, a + 774(0")) f) P(/ 3 ) = 0 . Hence ( 7 ) and ( 10 ) yield 

fcr+r 7 

11/ [/(y(t),t) - /(a;(t),t)]dt|| 
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< w(||W(cr+) — ^(cr)|| + a) 



P = P_ 

(n(/3) + l)w(||^(cr+) - ^(cr)|| + a ) n(/ 3) + 1 



( 12 ) 



for every 77 £ (a, a + 774 (a)). 

Denote U a (t) = n ^ +1 J2aeP(0) 1 <y( t ) for t G [to, ti] where / CT (t) = 0 for t < a and I a (t) = 1 
for t > a. The function U a : [to,ti] — > R is nondecreasing and continuous from the left and 



U a {t 1 ) - U a (t 0 ) = f n(p) < 0. (13) 

n((3) + 1 

The points of discontinuity of the function U a are clearly only the points belonging to P(l 3), 
and for t £ P(fi) we have 

U a {t+) — U a (t) = t £ [t 0 ,ti]. 

Using the function U ai we can set U a (t) = /3U c (t) + U a (t) for t £ [to,ti] where by U c denote 
the continuous part of the function U, then the function U a is nondecreasing and continuous 
from the left on [4,4] and by (9) and (13) we obtain 



U a (h) - U a (to ) = P(U c (ti) - U c {t 0 )) + U a (t 1 ) - U a (t 0 ) 

<P[U(t 1 )-U(t 0 ) + l) = ^. (14) 

If a € [ 4 , 4 ] \P(P) then set 6(a) = 773 (cr) > 0 and if a £ [to, 4] f) P(P) then set 6(a) = 

774 (a) > 0. By (11), (12) and definition of U a we obtain the inequality 

II [ [f(y(t),t)~ f(x(t),t)]dt\\<U a (a + r])-U a (a) 

J G 

for a £ [t Q , ti], 77 £ [0, , 5(a)], then by Theorem 1.3 and Theorem 1.17 in [5], we have 4>(|| f° +V [f(y(t), t) 
/(*(*), t)] ||) < &(U a (a + rj) - U a (a)) = 14(I/ a ; [a, a + 77 ]) < V$(U a ; [i 0 , cr + 77 ] ) - 14(14; [t 0 ,a]), 
therefor by (4) we have 



V(a + Ti,y(a + ii)) -V(a,x(a)) 



< V(V$(4'; [t 0 ,a + 77 ]) - 14(4'; [4, o']) + 14 (U a ; [ t 0 ,a + 77 ]) - 14 (U a \ [t 0 , a])) + r]R + 7 ?e 

= G(a+ 77 )-G(a), (15) 

for a £ [t 0 , 4 ], 77 £ [0,<5(a)] where G(t) = tV14(4'; [t 0 ,t]) + i?(t-to) + e(t-to)+tV14(/7 a ; [t 0 , t] ) , 
Gis of bounded <E>-variation on [to, ti] and continuous from the left on (4,4]. By (14), (15) and 
Lemma 2.2, we obtain 



4(4, 2/(4)) - V(t 0 ,y(t 0 )) < G(h) - G(t 0 ) 

= /V14(4'; [4, ti] ) + R(t\ — to) + s(t\ — to) + NV$(U a ; [to, ti]) 

< tV14(4'; [t 0 , ti]) + R(t\ — to) + e(t\ — 4) + N<5>(jj), 

Since e can be arbitrary small, we obtain the result given in (3) from this inequality. 

Theorem 3.2. Assume that V : [0, + 00 ) x B a — > R, 0 < a < c is such that for every 
x £ B a = {y £ R n : ||j/|| < a} the function V(-,x) is continuous from the left and 
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(1) Assume that the function V (t, x) is positive definite, i.e. there exists a continuous and 
increasing real function v : [0, +oo) — ► R Such that v(p) = 0 -<=>• p = 0; 

(2) for all (t,x) € [0, +oo) x B a , we have 

V{t,x) > t>($(M)) (16) 



and 

V(t, 0) = 0. 

For any x,y € B a and L > 0 is a constant, we have 



| V(t,x)-V{t,y)\ < L$(\\x - y\\). 



(17) 



(18) 



If the function V (t, x(t)) is nonincreasing and of bounded ^-variation along every solution 
x of (*) then the solution x = 0 of (*) is $-variationally stable. 

Proof. Since we assume that the function V(t,x(t)) is nonincrease and of bounded 
<F- variation whenever x : [cq/3] — > R n is a solution of (*), we have 



r V(t + y,x(t + ??))- V(t,x(t)) 

Iim sup 

77 ->- 0 + T] 



(19) 



for t £ [a, 0\. 

Let £ > 0 be given and let y : [to, ti] —l ► R n , 0 < to < t\ < +oo be of bounded <F-variation 
on [t 0 ,*i] and continuous from the left in (f 0 ,f i]- Since the function V satisfies the assumptions 
(2) of Theorem 3.1 with H = 0, We obtain by (3), (17) and (18) the inequality 



V(r,y(r)) < V(t 0 ,y(t 0 )) + NV^,((y(s) - J t S Q f(y(t),t)dt)-, [ t 0 ,r ]) 



< ^(||y(io)||) + NV^((y(s) - f t s g f(y(t),t)dt); [to, r-]) 

< N^(\\y(t 0 )\\)+NV^((y(s)-J t S o f(y(t),t)dt);[t 0 ,r}) (20) 

holds for every r £ [to, ti] - 

Define a(e) = inf r < e v(r), then a(e) > 0 for £ > 0 and lim e _ >0+ a(e) = 0. Further choose 
8(e) > 0 such that 2N8(e) < a(e). If function y satisfying 3>(||y(to)||) < 8(e) and V$((y(s) — 

/to /(y(i),t)dt); [t 0 ,r]) < S(e), then by (20), we have 

V(r,y(r)) < 2 N8(e) < a(e),r £ [t 0 , ti]- (21) 



If there exists a t £ [to, ti] such that < E>(||?/(t)||) > £ and by (16) we get the inequality 
V(t,y(t)) > u(<I>(||2/(f)||)) > inf r < e u(r) = a(e), which contradicts (21). Hence $(||j/(t)||) < £ 
for all t £ (to, ti] and by Definition 2.7 the solution x = 0 of (*) is $-variationally stable. 

Theorem 3.3. Let V : [0, +oo) x B a — > R be a function with the properties given in 
Theorem 3.2. If for every solution x : [to, ti] — > B a of bounded ^-variation of (*) the inequality 



lim sup 

t?— >0+ 



V(t + r], x(t + y)) - V (t, x(t)) 
V 



< —H(x(t)) 



( 22 ) 



holds for every t £ [to, ti], where H : R n — > R is continuous, H( 0) = 0 ,H(x) > 0 for x ^ 0, 
then the solution x = 0 of (*) is $-variationally-asymptotically stable. 
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Proof. From (22) it is clear that the function V(t,x(t)) is of bounded ^-variation and 
nonincreasing along every solution x(t) of (*) and therefore by Theorem 3.2, the solution x = 0 
is $-variationally stable. 

It remains to show that the solution x = 0 is 4>-variationally attracting. From the $- 
variational stability of the solution x = 0 of (*), there is a do G (0, $(a)) such that if y : 
— > R n is of bounded ^-variation on [t 0 ,ti] where 0 < t, 0 < t\ < +oo,y is continuous from 
the left on (t 0 ,t i] and such that $(||j/(to)||) < and V<s>((y(s) - J* q f(y(t),t)dt); [t 0 ,ti]) < d 0 , 
then 3>(||y(t)||) < $(a) for t G [to,ti], by the definition of $(u), we have ||y(t)| < a, i.e for 
every t € (to, t±] we obtain y(t) G B a . Let £ > 0 be arbitrary, from the ^-variational stability 
of the solution xEOwe obtain that there is a 5(e) > 0 such that for every y : [t 0 ,ti] — > R n 
of bounded ^-variation on [fo,ii] where 0 < to < t\ < +oo, y is continuous from the left on 
(*o, *i] and such that $(||y(*o)||) < 6(e) and Vq,((y(s) - / t ® f(y(t),t)dt ); [t 0 ,t i]) < d(e), we have 
^(lls/WII) < £ f° r ^ € [to,£i]- Define 7 (e) = min(<5 0 , d(e)) and T(e) = _ j\r ' 5 o+ g'l e ) > 0 where 
R = sup {—H(x) : 7 (e) < $(||x||) < e} = — inf {H(x) : 7 (e) < $(||x||) < e} < 0 and assume 
that y : [to, ti] -> R n of bounded ^-variation on [t 0 , ti] where 0 < to < ti < + 00 , y is continuous 
from the left on (t 0 , ti] and such that < h(||j/(to)||) < S 0 and 

V<s>((y(s)~ f f(y(t), t)dt); [t 0 , ti]) < 7 (e). (23) 

Jto 

Assume that T(e) <t\—to,i.e.to+T(e) < t\. We show that there exists a t* G [to, to+T(e)] 
such that d>(||y(t*)||) < 7 (e). Assume the contrary, i.e. 4>(||y(s)||) > 7 (e) for every s G [to, to + 
T(e)]. Theorem 3.1 yields 

V(t 0 + T(e),y(t 0 + T(e))) — V(to,y(to)) 

< NV&((y(s) - f(y(t),t)dt ); [t 0 , t 0 + T(e)}) + RT(e) 

< N"/(e) + = -N5 0 . 

Henc eV(to+T(e),y(to+T(e))) < V(to,y(t 0 ))-N6 0 < L$(||2/(to)||)-^o < N$(||y(t 0 )||)- 
N6 0 < (N — N)S 0 = 0, and this contradicts the inequality V(t. 0 + T(e),y(t 0 + T(e))) > 
u(<I>(||y(to + T(e))||)) > v(j(e)) > 0. Hence necessarily there is a t* G [to, to + T(e)] such 
that $(||y(t*)||) < 7 (e). And by (23), we have 3>(||y(i)||) < e for t G [t*,ti]. Consequently, 
$(||y(t)||) < e for t > t 0 + T(e). Therefore, the solution x = 0 is a $-variationally attracting 
solution of (*). 

Remark 3.1. If the function <f>(w) is defined in last section such that 0 < < +00 then 

by Theorem 1.15 in paper [5], we have BV<s>[a,0\ = BV[a,0\ where BV$,[a,/3] and BV[a,/3\ 
denote the classes of the functions bounded ^-variation and the functions bounded variation 
on [a, /3] respectively in usual sense. Hence, the main results of this paper are equivalent to the 
results of variational stability in [4]. 

If lim tI _ >0 + = 0 then by Theorem 1.15 in [5], we have BV[a,(3\ C BV§[a,0\. Such as if 
$(u) = u p ( 1 < p < + 00 ) we have lim u _ >0 + ^ 77 ^ = 7 T = 0- Therefore these results are essential 
generalization of variation stability of bounded variation solution to discontinuous system in 

[4]- 
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§1. Introduction 

The Pell and Pell-Lucas sequence satisfy the following recurrence relations, respectively: 

Pn = 2P n _i + P n - 2 , P 0 = 0,P\ = 1; 

Qn — ^Qn—1 T Qn—2>Qo Ql 2. 

Their n th are given by 

n - P n 

2V2 ’ 

Q n = a n + 0 n . 

where a = 1 + \/2 and 0 = 1 — \/2. 

Bozkurt I 1 ! constructed the right circulant matrices with Pell and Pell-Lucas sequences. 
These matrices take the following froms, resepectively: 





< Po 


Pi 


P 2 • 


Pn-2 


Pn- 1 \ 




Pn- 1 


Po 


Pi ■ 


Pn-3 


Pn-2 


RCIRC n (P) = 


Pn- 2 


Pn- 1 


Po ■ 


Pn- 4 


Pn-3 




P 2 


P 3 


P 4 . 


• • Po 


Pi 



V Pi P 2 P 3 -. Pn - 1 P 0 / 



( 1 ) 

( 2 ) 
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RCIRC n {Q) = 



Qo 


Qi 


Q 2 


Qn—2 


Qn—l 


Qn—l 


Qo 


Qi 


Qn — 3 


Qn—2 


Qn — 2 


Qn—l 


Qo 


Qn — 4 


Qn - 3 


Q 2 


Q 3 


Qi 


Qo 


Ql 


Ql 


Q 2 


Q 3 


Qn—l 


Qo 



He presented the determinants and the ineverses of these matrices. 



( 4 ) 



§2. Preliminary results 



The following will be used to prove some of the main results. 

Lemma 2.1. 



n— 1 

£ [«* - />*] “ 
k—0 

where u> = e 2m ^ n . 

Proof. 

n— 1 

£[“ fc -0 fc ] w 

k = 0 



—mk 



a n — /3" + (a n ~ 1 - p n -Quj- m -(a-p) 

w -2m + 2u}~ m - 1 



—mk 



1 - a r ‘ 



1-/3" 



1 - aor m 1 - pu)~ m 

(1 - a n )(l - /3 Lo~ m ) - (1 - /3 n )(l - au~ m ) 

1 - 2w~ m - oj~ 2m 

a n - p n + (a"" 1 - /3"- 1 )u;- ro -{a-p) 

w -2m + 2 w - m - 1 



Lemma 2.2. 



n— 1 



53 [a fe - /3 fe ] 2 = Q2n + g2 "~ 2 - (-1)". 



fc— 0 



Proof. 



n— 1 



n— 1 



53 [a fe - ^ fc ] " = 53 [ a2k - 2 ^) k + P 2k 



k — 0 



k = 0 



1 — a 2n , l-^ 2n _ 



1 - a 2 1 - /3 2 “V 2 
(1 - a 2 ")(l - /3 2 ) + (1 - /3 2n )(l - a 2 ) 
(1 — a 2 )(l — P 2 ) 

-4 - (a 2n + p 2n ) - (a 2 ”" 2 + ^ 2 ”- 2 ) 



+ (-l) n -l 



Q2n + Q2n-2 



-4 

-(-ir 



+ (-!)” -1 



Lemma 2.3. 



n— 1 



53 [ a fc + ^ 2 = Q2ra + Q2n - 2 + 2 - (-1)". 



k—0 



( 5 ) 



( 6 ) 



( 7 ) 
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Proof. 

n— 1 

E [“*+£*] 2 = 

k — 0 



For the rest of the paper, we will use |A|, ||^4 ||b and ||A ||2 to denote the determinant, 
euiclidean norm and spectral norm of matrix A , respectively. 



n— 1 

E [« 2fc + 2(a/3) fe + /3 2 1 

fc= o 

1 - a 2n 1 - /3 2 " 



1 — a 2 






+ 2 






(1 - cr")(l - /3 2 ) + (1 - /3 2n )(l - a 2 ) 
(l-a 2 )(l-/? 2 ) 

-4 - (a 2n + ^2”) - (a 2 "" 2 + ^ 2n - 2 ) 



+ l-(-l) r 



Q 2n + Q2n-2 



-4 

+ 2 — (— 1)". 



+ i-(-ir 



§3. Main results 



Theorem 3.1. The eigenvalues of RC I RC n (P) are given by 



Proof. 



^ m — 



(Pn - 1 - 1 )cQ- m + Pn 

UJ-2 m + 2ui~ m - 1 



n— 1 



^ ^ Pk W 



— mk 



k—0 



= E 



[a fc - P k }io- mk 

2^2 



k = 0 



a n - /3 n + (a"" 1 - f3 n - 1 )u~ rn -(a-/3) 

2 v / 2(w- 2m + 2w m - 1) 

(Pn-! - 1 + Pn 

OJ-2 m + 2uj~ m - 1 



( 8 ) 



Theorem 3.2. The eigenvalues of RC I RC n (Q) are given by 



Qn+i + 2(cc- m - 1) 
w -2m + 2u>~ m - 1 



( 9 ) 
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Proof. 



Theorem 3 



Proof. 



Theorem 3 



Proof. 



Theorem 3 



n— 1 

Hm = 

k—0 
n— 1 

k = 0 

1 - a n 1 - /3 n 
= 1 — 

1 - au- m 1 - /3w" m 

(1 - a n )(l - /3 uj~ m ) + (1 - /3 n )(l - au;-™) 

1 - 2u}~ m - oj~ 2m 

2 -Q n - 2u>~ m - 

1 - 2uj~ m - UJ~ 2m 

Qn + 1 + 2{u~ m - 1 ) 

w -2m + 2u)~ m - 1 



RCIRC n (P) 



\j2n( y Q-2n + Q 2 n- 2 ) ~ (— l) n 64 
8 



RCIRC n (P) 



n— 1 




a fe — f3 k 

2y/2 




P k } 2 



J n 


Q2n Q2n—2 ( 1 




Vs 


4 ^ c 




yj2n(Q2n + Q^n— 2 ) — ( — l) n 64 



8 



4. 



\/ n(Q 2 n + (? 2 n— 2 ) + 8 — (— 1)"4 
2 



RCIRCn(Q) 



n— 1 



[a fe +/3 fe ]" 

\ fe=0 



= \ n 



Q2n + Q 2n-2 



+ 2 — (— l) r 



\/ n(Q2n + Q2n-2) + 8 — (— 1)"4 



.5. 



■Fn-l ~ Pn ~ 1 



< 



RCIRC n {P) 



< 



Pn+l ~ 1 



(10) 



( 11 ) 



( 12 ) 




Vol. 10 



Note on right circulant matrices with Pell and Pell-Lucas sequences 



15 



Proof. Note that 



RCI RC n (P) = max { j A m | } , 



hence, we have the following inequality 

|(-Pn— 1 - l)w" m | -P, 



w 



+ 2oj~ m - 1 | 

Pn - 1 — Pn — 1 



Pn-l - Pn ^ 1 



< 



< 



< 



RCIRC n (P ) 
RCIRC n (P) 
RCIRC n (P) 



< 



|(P»-! ~l)iQ- m \+P n 

\ UJ -2m _|_ 2uj~ m - 11 



< 



< 



Pn-l — 1 P Pn 



Pn+1 — 1 



Theorem 3.6. 



Qn—l 



< 



RCIRC n (Q) 



< 



Qn+l + 4 



Proof. Similar as the previous theorem 3.5. 



(13) 



References 

[1] D. Bozkurt, On the determinants and inverses of circulant matrices with Pell and Pell- 
Lucas numbers, arXiv: 1201, 6061vl [math. NA], 29 Jan, 2012. 




Scientia Magna 
Vol. 10 (2014), No. 1, 16-33 



Evaluation of a summation formula confederated 
with contiguous relation 

Salahuddin 
Mewar University, 

Gangrar, Chittorgarh (Rajasthan) , India 
E-mail: vsludn@gmail.com sludn@yahoo.com 

Abstract The main aim of the present paper is to evaluate a summation formula confederated 
with contiguous relation and recurrence relation. 

Keywords Gaussian hypergeometric function, contiguous function, recurrence relation, Bailey 
summation theorem and Legendre duplication formula. 

2010 Mathematics Subject Classification: 33C60, 33C70, 33D15, 33D50, 33D60. 



§1. Introduction and preliminaries 



by 



or 



Definition 1.1. Generalized Gaussian hypergeometric function of one variable is defined 



aFb 



ai, a.2, ■ ■ ■ ,a,A 



z 



bi,b2, • ■ ■ ,bs ; 




(«l)fc(Q2)fc • • • ( dA)kZ k 

(bi)k(b 2 )k • • • (^s)fcfc! 





Wa) ; 




1 

5! 

1 


aFb 


Z 


= aFb 


Z 




(b b ) ; 




1 

oql, 

1 




{{a A ))kZ k 



( 1 ) 



where the parameters b\,b 2 ,- ,6b are neither zero nor negative integers and A , B are non- 

negative integers. 

Definition 1.2. Contiguous relation [E. D. p. 51(10), Andrews p. 363(9. 16)] is defined 
as follows: 



(a - b ) 2 Ei 



a, b ; 




a + 1, b ; 




a, b + 1 5 


z 


= a 2 Ei 


z 


1 

O- 

to 


z 


c ; 




c ; 




c ; 



Definition 1.3. Recurrence relation of gamma function is defined as follows: 



( 2 ) 



T(z + l)=zT(z). 



( 3 ) 
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Definition 1.4. Legendre duplication formula [ Bells & Wong p. 26(2.3.1) ] is defined 
as follows: 



y/n T{2z) = 2 (2z ~ 1) r{z)r(z+^j, 

2 (b-i) r (6) r (6±l) 



r ( - 1 = Vtt = 



m 



( 4 ) 

( 5 ) 

( 6 ) 



_ 2 (°- 1 ) r(§) r(s±i) 

Definition 1.5. Bailey summation theorem [Prud, p. 491(7.3.7.8)] is defined as follows: 



2-F1 



a, 1 — a ; 1 
• 2 



r(|) r(^) 

r(c±a) r(£±iz^) 2'- 1 r(!f)r(^)' 



T(c) 



§2. Main results of summation formulae 



( 7 ) 



2 -Pi 



a , — a — 49 ; 1 

. 2 



/n r(c 

2 c+49 



40015948329812641136799430213632000000 



P(C^) r (£+^f49) 



-+ 



+ 



2 / \ 2 

-74697042588755634816870901715435520000a 

r(^) r(^t®j H 

45701551565700312458589664488136704000a 2 

p(c-a) r (c±a±49) 4 

-12557326585358016296749495447708262400a 3 

P(c-a) r (c+<|t49) 

| 1607017301962727354472466577244395520a 4 i 

+ r(^) r(^®) + 

-66334092838839671097859298211598848a 5 - 4173963485546119223032377440607360a 6 

r(s=2) r(£±^t«) 

354268334898336692529886398612096a 7 + 6493428631831030022245150741472a 8 



+ 



- + 



+ - 



P(c^a) r(c+^t49) 

-716075598615688398506581518048a 9 - 13409906223325457338650846040a 10 

r(^) r(^t«j 



-+ 



695228979926866432431858936a 



11 



20318017073982774702824162a 



12 



P(c^a) r (c±^f49) 

-181081529755913127436368a 13 - 12758198887309415982820a 14 



P(c^a) r(£±g_H9) 
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— 113809757878112295384a 15 + 1948115287787037182a 16 + 46874133673870272a 17 

r(£=«) r(£±s±M) 

257645168805200a 18 - 1896116954424a 19 - 31420958338a 20 - 136837008a 21 

J L 



+ 



r( 



c—a 

2 



) r( c+a + 49 ) 



-32980a 22 + 1176a 23 + 2a 24 + 102873723702221439014473177529057280000c 

r(s=“) r(£±^9) 

— 143953054337702235151044038042320896000ac 

+ r(^) r(£±^l9) + 

69785186656677903329815708782664089600a 2 c 



+ 



+ 



r(c^a) c+(^f49^ 

— 152540878242863373435008530 16781 127680a 3 c 



T(c^a) r(£+£Lt49) 



; 1470184367841770478008747704464310272a 4 c ; 

+ r(^ ) r(£±^) + 

— 25186819706236967969752199080869888a 5 c — 4916413813362186868890467613499392a e c 

r(s=2) r(£±^) 

180523758084624041877205846241280a 7 c+9444650003452625898653645586432a s c 

J 



+ 



r( 



c—a 

2 



) r( c+a + 49 ) 



+ 



— 290073752007187728910991044608a 9 c — 13790932677896490825858473472a 10 c 
+ r(s^) r(s±^t49) 

i 132168254665462941656989440a 11 c + 11478210451803468174127872a 12 c ( 

+ r(^) r(£±^H) + 

82419396059301592042752a 13 c - 3441243829617943295232a 14 c 



r( 



c —a 

2 



) r( c+a + 49 ) 



— 68832268710004800000a 15 c - 120844174461424128a 16 c + 874343974491 1872a 17 c 

re n r ■ 'v ;,, i 

; 96443587688448a ls c + 222439176960a 19 c — 2314760448a 2O c — 16144128a 21 c ; 

+ r(^) r(£±^l9) + 

— 29952a 22 c + 107759352272213925707820568758190080000c 2 



r( 



c—a 

2 



) r( c+a + 49 ) 



+ 



+ 



— 119604791694733248989672753097867264000ac 5 
r(s=“) r(£±^9) 

; 46915581673954062917864644291866132480a 2 c 2 

r(^) r(£±^M) 

— 8179716567466648414966719773717397504a 3 c 2 



+ 



r( C -T) r( 



c+a+49 



| 572414151475230486442581108591747072a 4 c 2 i 
+ r(^) r(£4^ti9) + 
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4225122815762079559693237812774912a 5 c 2 - 2107743828904816275105349039067648a e c 2 

21777324127000959865037942782208a T c 2 + 3811313237604914006353415570432a 8 c 2 

+ r(£fi) r(£±^) 

- 17208144178264485968744984832a 9 c 2 - 4072656512763355130085794976a lo c 2 



-+ 



r(^) r( 



c+a+49 



— 30448585606100284800091824a 11 (r + 1971285851965809130849552a 1 V 



W) r( 



c+a+49 



+ - 



37909392941438882335472a 13 c 2 - 136814112867707752048a 14 c 2 



r(c-a) r ^c+a^f49^ 

— 100 18332037253860320a 4 5 c 2 - 86695521821402208a 16 c 2 + 228479601333984a 17 c 2 



r(c^a) r (c+^±49) 

7725554252416a 18 c 2 + 40796996240a 19 c 2 + 25577552a 20 c 2 - 336336a 21 c 2 



- + 



T(c-a) r (c+i+t49) 

— 624a 22 c 2 + 64423398177963803144998543335358464000c 3 

r(£^) r(£±^) 



-+ 



+- 



-58300150202899337016565241575133675520ac 3 



T(c-a) r (c+<|t49) 



18691960982318475085206332143677997056a c 



2„3 



r(c=a) r (c+^±49) 



+ - 



— 2586065877925151796735934133731590144a 8 c 3 



' c+a+49 



2 , 

122171345804374971207330799098003456a 4 c 3 



1 r(^) r(s±^) + 

4533540855919510991454734335180800a 5 c 3 — 467631969431579769690069597929472a e c 3 



+ - 



r(c^a) r (c+^±49) 

-4697327171517980882726705283072a T c 3 + 722969542181022180843132002304a s c 3 



+ - 



T(c^a) r (£i+ti9) 

9083202730777950404000716800a 9 c 3 - 52957857198852051 1523748864a 10 c 3 

r(^) r(£±^) 



— 11402466632175022106486784a ii c 3 + 108169499556766290779136a 12 c 3 



r(c=a) r (c+^±49) 



- + 



5063222588621911633920a i3 c 3 + 32000805092958873600a 14 c 3 - 498995580905902080a 15 c ; 



,15„3 



W) r( 



c+a+49 



— 8205717988435968a 16 c 3 - 29314507622400a 17 c 3 + 167330042880a ls c 3 + 1506785280a 19 c 3 

F( ' ) lVV 19 ) 

3075072a 2O c 3 + 25347138973058359260380026128826368000c 4 
+ r(£=“) r(£±#3lj + 
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Salahuddin 



No. 1 



+ - 



-18978543797716490715225178573727858688ac 4 

r(^) r(£±^H) 

4991684025479138391935744923474722816a 2 c 4 



+ - 



r(c-a) r (c+£oH9) 

— 541099749641591095487701589399519232a 3 c 4 



r(c=a) r( cM49) 



+ - 



14606689411535623104290723942940672a c‘ 



4 „4 



' c±a + 49 



ww 2 , 

1292549598384288914173152912040960a 5 c 4 - 59385691383614919077242368788992a e c 4 

IVY'! IV "V ; "i 



—1892451 113166353101369696786944a'c‘ 



7^4 



72169632318929524275249528576a 8 c 4 



F(c=a) r (c+^49) 

2228687725107369556056237120a 9 c 4 - 28472710924167755430485664a 10 c 4 



_r r(s=2) r(£±^t«) + 

— 1409430360340354713769408a 11 c 4 - 5501770617352422372960a 12 c 4 

+ r(£=“) r(£±^) + 

297674907085096661 120a 13 c 4 + 3995332631611006400a 14 c 4 + 1153916414851200a 15 c 4 

r(^) r(£±^9) 

— 295200075538368a 16 c 4 - 1932687556800a 17 c 4 - 1947065120a 18 c 4 + 15695680a 19 c 4 
f r(^) r(5±^ 49 ) 

32032a 2O c 4 + 7099535118854336323202874869003845632c 5 

+ r(^) r(£±^l9) + 

— 44323362301 10420896646668265044574208ac 5 



- + 



F(c^a) r(c+^t49) 



2 / V 2 

955466413384038667615798431427067904a 2 c 5 



- + 



r( c=a) r ^ ^ 

— 7938039714 1399923587473758714593280a 3 c 5 -t - 64380593 162 16642848934951 24295680a 4 c 5 

r(£^) r(£±^l9] 

209465328810891152290757390991360a 5 c 5 - 4028319000276505780447170035712a 6 c 5 



+ - 



+ - 



r(c-a) r ^c+(|f49 j 

— 2909 15966729411 309455672688640a 7 c 5 4- 3 124145533 11 13666 14800302080a 8 c 5 



2 > ' \ 2 

241803413658639276317153280a 9 c 5 + 584445378775657625081856a lo c s 



W) r( 



c+a+49 



- + 



— 87844798556732071895040a i V - 1067805772051345367040a 12 c 5 



F(c-a) r (c+a±49) 



- + 



6128657556733071360a 13 c 5 + 189548697075978240a 14 c 5 + 958658985123840a 15 c 5 

r(^ ) r(£±#^9) 
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— 3202625986560a ie c 5 - 40683202560a 17 c 5 - 92252160a 18 c 5 

+ r(^) r(£±^49^ + 

1487119300880904242605217024921567232c 6 - 776777553617389082534820001691467776ac 6 

r(£=“) r(^"- 40 ) 

136593392051968205997829269713846272a 2 c 6 - 8388092183067964659019677414162432a 3 c 6 

r(^) r(5±^t«j 

-91340670507632683968759660462080a 4 c e + 22224973282782547275305692944384a 5 c 6 



- + 



' c±a + 49 



+ - 



r(^)ro 2 , 

— 35529674554603169040172140544a e c 6 - 26210702433049018598332609536a 7 c 6 

r(c-a) T (c+a_ t49) 

- 103200723777727230528706560a 8 c 6 + 14895636537534753758389632a 9 c e 



+ - 



F(c=a) r (W49 ) 



'c+a+49'i 
>. 2 > 

182403094584283413160320a 10 c 6 - 2663633385875852662272a n c 6 
+ r(£=fi) r(s±^Mj + 

— 62256983715951664640a 12 c e - 174817816321564416a 13 c e + 4084756093786368a 14 c e 

r(s=2) r(£±^9j 

34430420344320a 15 c 6 + 48609200640a 16 c 6 - 282522240a 17 c 6 - 640640a 18 c 6 



+ - 



r(c^a) r(c+^t49) 



240838218499636521398013613710508032c 7 - 105287696394810237973591256098406400ac 7 



r(c=a) r(^ 49 ) 

14976654351475036084045130810523648a 2 c 7 - 638818690035905266292950472589312a 3 c 7 

r(c-a) r (c+^f49) 

— 20960541592500617722444448268288a 4 c 7 + 1615109961509085989321865953280a 5 c 7 | 
+ r(^) r(£4^ti9) + 

21023698439479542356321501184a 6 c 7 - 1497276217131824539913158656a 7 c 7 



+ 



r(c^a) r( c±<^49 ) 



- + 



+ - 



— 23179379774042918881787904a 8 c 7 



515692831831999968706560a 9 c 7 



r(c^a) r(c±^t49) 



12163046240440205180928a lu c 7 - 6024819137364099072a 14 c 7 - 1837489800964866048a 12 c 7 

r(^) r(£±^i^ 4 

— 12946877731307520a 13 c 7 + 23007583272960a 14 c 7 + 495947612160a 15 c 7 + 1265172480a 16 c 7 

r(^) r(5±^t®j 

30871602486909275910830889875013632c 8 - 11272571429148914722533062085705728ac 8 



' c±a + 49 



WW 2 , 

1280968945576182794260760676794368a 2 c 8 - 33589734750958810846109982916608a 3 c s 

iv ./'! rr"V ; "i ' ‘ 

i — 2206717841266713417081546031104a 4 c 8 + 79728366271118889239432724480a 5 c 8 | 

+ r(£=“) r(£±#^9) + 




22 



Salahuddin 



No. 1 



2306069297890733863547265024a 6 c 8 - 52263325577453158896279552a 7 c 8 

+ r(^) r(£±^49^ + 

— 1561668419447767019271168a 8 c 8 + 6335745740330385899520a 9 c 8 
+ r(£=“) r(£4^ti9) + 

433099863114577072128a lo c 8 + 2445053402814038016a n c 8 - 26209896905127936a 12 c s 

r(^) r(5+^®) 

— 303502718238720a 13 c 8 - 568422666240a 14 c 8 + 2583060480a 15 c 8 + 6589440a 16 c 8 



' c±a + 49 



WW 2 , 

3185642085289031342835468884508672c 9 - 967550638958541171077959970193408ac 9 
’ " l'C ./'! lev 

86299774734978829551171583082496a 2 c 9 - 998945521658261204553034629120a 3 c 9 
+ r(^) r(^±^49) + 

-153617690085710615493130321920a 4 c 9 + 2411197000642734061356318720a 5 c 9 



+ - 



r(c^a) r (£±^f49) 

137452472767000072564506624a e c 9 - 785064641247377381 130240a 7 c 9 

+ r(£fi) r(£±^) 

— 61009257199589456609280a s c 9 - 256032809547607572480a 9 c 9 

+ r(^) r(£±^i9j + 

8769591 108516446208a lo c 9 + 88785671187333120a u c 9 - 46103447470080a 12 c 9 



r(c=a) r (c±^±49) 



- + 



+ - 



— 3278764769280a 13 c 9 - 9559080960a 14 c 9 + 267889258438085483657974117302272c 10 

r(c^a) r ^c+^ti9) ^ 

-67263084939581397145810722881536ac lo + 4593733502680143944975661400064a 2 c 10 

r(£=“) r(£±^9) 

8539339140610043749157830656a 3 c 10 - 7676662187554452553435054080a 4 c 10 



F(c=a) r( W49 ) 



^ c+a+49 \ 

>. 2 > 

19441815522855857897013248a 5 c 10 + 5387422458233145326796800a 6 c 10 

r(£ft) r(£±^H) 

20474586985880613634048a 7 c 10 - 1490183150114131066880a 8 c 10 



-+ 



T(c-a) r (£+<|t49) 



- + 



— 145432505261 16888576a 9 c 10 + 83327898738180096a 10 c 10 + 1503081400877056a n c 10 

r(^) r(5±^®j 



-+ 



3662123745280a 12 c 10 - 13661519872a 13 c 10 - 39829504a 14 c 10 



' c±a + 49 



- + 



r (^W 2 , 

18518169236138399504184033411072c 11 - 3812360305782072606527229788160ac n 

F(c-a) r (c+a±49) 



1923406886849702093574541 14816a 2 c 



2„11 



2855197808016937994321657856a 3 c 11 



r(c^a) r(£±^) 



-+ 
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+ - 



— 281873037311025074002722816a 4 c 11 - 2077996489024813179863040a 5 c 11 

r(£=2) r(^±^9j + 

i 143506954311464043675648a e c 11 + 1541710276583120437248a 7 c n | 

+ r(^) r(£4^ti9) + 

-21366367851509710848a 8 c n - 341036842572840960a 9 c n - 235238901940224a 10 c n 



r(c=a) r( cM«) 



+ - 



12876603457536a 11 c 1 



43797970944a 12 c 1 



1058340146221982392995043868672c 1 



r(c-a) r ^c+^ 4 y^ 

-176775637161072795796969095168ac 12 +6235849644708280985967394816a 2 c 12 

re./'! ic v ; "i 



+- 



179112003462820022860382208a 3 c 12 - 7549631558688412740485120a 4 c 12 

r(^) r(5+^ 49 ) 

— 121866570019913213214720a 5 c 12 + 2500009663439970992128a 6 c 12 



+ - 



+ r(£=2) r(^t49) 

44613295101771841536a 7 c 12 - 113644325380325376a 8 c 12 - 4446228579287040a 9 c 12 

r(^ ) r(£±^) 



— 14131384909824a 10 c 1 2 + 44710428672a n c 12 + 152076288a 12 c 1 
r(£=“) r(£±^9) 



-+ 



50169135355853067670825992192c 13 - 6709363518959579397365956608ac 13 



r(c=a) r (c±^±49) 



- + 



+ - 



150614391320542327439622144a 2 c 



2„13 



.13 



F(c-a) r (c±<^4y) 

— 140953864332077212631040a 4 c 13 - 3711753500792413224960a 5 c 13 



+ 



+ r(^) r(£±^t«) 

23530276519912931328a 6 c 13 + 762927749215027200a 7 c 13 + 1587081243525120a 8 c 13 

rx^mra 9 ) 

— 31366085345280a 9 c 13 - 128024838144a 10 c 13 + 1974250608891543017637281792c 14 

r(^) r(£±^H) 

— 207912674105719908850466816ac 14 + 2433971414854765171965952a 2 c 14 



+ - 



r(c-a) r (c+<|t 49 ) 

197419671000565497200640a 3 c 14 - 1569222994723146301440a 4 c 14 



r(c=a) r (c+<^49 ) 

— 73662625200783360000a 5 c 14 - 28046683124662272a 6 c 14 + 8006131050086400a 7 c 14 



W) r( 



c+a +49 



33986276229120a 8 c 14 - 93351444480a 9 c 14 - 381026304a lo c 14 ; 

+ r(^) r(£±^) + 

64395185753779071029870592c 15 - 5229548355031289718374400ac 15 
+ r(^) r(£4^ti9) + 




24 



Salahuddin 



No. 1 



14788626650875350220800a 2 c 15 + 4080138857445379276800a 3 c 15 

+ r(^) r(£±^49^ + 

— 1838738814365859840a 4 c 15 - 984309579582013440a s c 15 - 3963648988938240a e c 15 

r(s=2) r(s±^j 

47795939573760a ? c 15 + 243856834560a s c 15 + 1733754269178556886024192c 16 



r(c^a) r( c±<^49 ) 

— 105699998215764315209728ac 16 - 462776974461593190400a 2 c 16 



+ 



F(c^a) r(c±^ t49) 

62047273223138574336a 3 c 16 + 274014892392579072a 4 c 16 - 8585594583121920a 6 c 16 

T(c-a) T (c+a_ H9) 



+ - 



-51289950781440a 6 c 16 + 124468592640a 7 c 16 + 635043840a 8 c 16 



f c-j-a+49 



n^)r^ 2 , 

38253074247649509507072c 17 - 1690248262256299081728ac 17 



^ c+a+49 



— 17030905158269140992a c 



2J7 



r (^) T (~ 2 7 

677374598861291520a 3 c 17 + 5103823135703040a 4 c 17 



r(c=a) r (W49 ) 



' c+a+49 ^ 
v 2 > 

— 44281532252160a 5 c 17 - 301234913280a 6 c 17 + 683995510182533660672c 18 



r(c=a) r(? M49) 



’ c+a+49 'i 
>. 2 > 

— 20883614573480902656ac ls - 300907292140765184a 2 c 18 + 5031796612792320a s c 18 



+ - 



F(c-a) r (c+<|t49) 

47159301898240a 4 c 18 - 102503546880a 5 c 18 - 697303040a 6 c 18 + 9749520870090473472c 10 

r(s=2) r(£±s±®) 

— 192157335441899520ac 19 - 3363268429086720a 2 c 19 + 22788156948480a 3 c 19 

+ r(£ft) r(£±^l9) + 

232532213760a 4 c 19 + 108101378949251072c 20 - 1238931957874688ac 20 

+ r(^) r(£±s±M} + 

— 24121 180160000a 2 c 2 ° + 47475326976a 3 c 20 + 484442112a 4 c 20 + 898096988946432c 21 | 

r(5=«) r(e±^ 49 ) + 

— 4991691522048ac 21 - 101871255552a 2 c 21 + 5257174188032c 22 - 9453961216ac 22 



r(c^a) r(c+^f49) 

— 192937984a 2 c 22 + 19327352832c 23 + 33554432c 24 

F(c-a) T (c+a_ H9) 

-174640699428026688529844729138380800000a 



r (c_a+l) r( - c+a+50 ^ 

203146678984534054497568403939819520000a 2 



r (c_a+l) r ( c+a+50 ) 

-85138755175820398653100504242008064000a 3 

r(g=fS) r( c+a + 50 ) 



-+ 
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+ - 



16061280659531143673149881316338278400a 4 



r (c_a+l) r ( c+a+5U ) 

—1222084280207838423985096126241433600a 5 



r (c_a+l) F ( c+a+50 ) 

-12283691196452377116468673849996800a 6 + 5784460407216745090330911562252800a 7 



r (c_a+l) r ( c+a+5U ) 

-64163356321096906750752842736000a 8 - 13694126502362695192734160137600a 9 



- + 



F (c_a+ 1 ) F( . c+a+50 ) 



70187206573178654636297887200a 



10 



19908972323990635614823432800a 1 



r (c_a+i) r ( c+a+5U ^ 

173269159811839104008127000a 12 - 13877737266129829896585600a 13 

r(g=|+i) r( c+a + 50 ) 

-321996506691917154610800a 14 + 1564939330380426184800a 15 



+ - 



137080053389848041000a 16 



r( c+a + 50 ) ~ + 

1531503700614518400a 17 - 5877737170852800a 18 



r( 



c— Q+l 



) r( 



c+a+50 



2 / V 2 / 

-269958340159200a 19 - 2165643711000a 20 - 2271561600a 21 + 62053200a 22 

r(g=|+i) r( c+a + 50 ) 

352800a 23 + 600a 24 + 174640699428027308978246462377820160000c 

+ r( c+a 2 f50 ) + 

— 532864667571457665192973079183032320000ac 

+ r(s^|+i) r( c+a + 50 ) + 

402516691346265375314108597722509312000a 2 c 



F (c_a+ 1 ) F ( c+a+50 ) 



- + 



— 125889625791092176675310990146311475200a 3 c 



+ - 



F (c_a+ 1 ) F ( c+a+50 ) 

17789322238914085698760878421723806720a 4 c 

r(^|+i) r( c+a + 50 ) 

— 839642652359017269802420497774739200a s c 



- + 



r(^=f±i) r( c+a + 50 ) 



-+ 



—437377996641 11358026071194502939200a 6 c + 4327054136964868650145069269816000a 7 c 



+ - 



F (c_a+ 1 ) F ( c+a+50 ) 

657101 16059613413534021275073520a 8 c — 8866517878966272728311508397200a 9 c 



- + 



F (c_a+ 1 ) F ( c+a+50 ) 

-150823836484592478601551835500a lo c + 8809059086367927385757986500a n c 



F (c = a+l) F ( c+a+50 ) 

245730574725150062467057345a 12 c - 2419735295216118363229200a 13 c 



r( 



c— a+1 



) r ( 



c+a+50 



-+ 
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— 158884196614329578449050a 14 c - 1380872827191951285500a 15 c 

+ r(^|+i) r( c+a + 50 ) + 

24710949541561727895a 16 c + 584963403786946800a 17 c + 3188363454621400a ls c 



— 23871357856500a 19 c - 392868365505a 20 c - 1709061200a 21 c - 409650a 22 c 

r( c^i+i) r( c + a + 50 ) 

14700a 23 c + 25a 24 c + 329717988586925953482621073961779200000c 2 



—5936101 15863251160379838852807786496000ac 2 



332339392507645324571342448478991155200a 2 c 2 

r( c=|+i) r( c + a + 50 ) 

-80292001276403892913036261386947788800a 3 c 2 



8435752253978765649343405490604441600a 4 c 2 



-182800700903842751424899605150310400a 5 c 2 



-27864925923991726587974486723020800a 6 c 2 



— 1866238092614775212796148761600a 9 c 2 - 82875721352461554443602089600a lo c 2 



896036369449822368051211200a n c 2 +71047302433059217158129600a 12 c 2 

+ r(g=|+i) r( c+a + 50 ) + 

4835662950387921 18753600a 13 c 2 - 21705016205569712558400a 14 c 2 

+ r(g=|+i) r( c+a + 50 ) + 

— 427069008643285008000a 15 c 2 - 704037139627536000a 16 c 2 + 54829120781635200a 17 c 2 



601725264940800a 18 c 2 + 1380827448000a 19 c 2 - 14486472000a 29 c 2 - 100900800a 21 c 2 



— 187200a 22 c 2 + 276232179692810109214204037425004544000c 3 



— 361016488745751900835386008702130585600ac 3 



157654746786562756857346716177587896320a 2 c 3 
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—2981 1845754423340406338901058806169600a 3 c 3 



r (c_a+l) r( . c+a +5U ) 



- + 



+ - 



2261040469101812046662593980909465600a 4 c 



4 „3 



F(c^a+1) 



c+a+50 



2 / \ 2 / 
7467962566571130491482924613606400a 5 c 3 



r (c_a+l) F ( c+a+50 ) 

-8288810131500283801973700171383040a 6 c 3 + 105753830113991285004211505699200a 7 c 3 



r( 



C— Q+l 



)T( 



c+a+50 



-+ 



+- 



2 7 V 2 1 

15276314790358723950186005280000a s c 3 — 89839656535308780981746947200a 9 c 3 



r( 



c— a+1 



) r ( 



c+a+50 > 



2 / V 2 / 

— 16726883912057639605189273680a lo c 3 - 115201961703676996712982600a n c 3 



r (c_a+l) r( - c+a+50 ) 

,12 + 



+ - 



8257727829380832345907000a 12 c 3 + 155470133092201748020200a 13 c 3 



F (c_a+ 1 ) F ( c+a+50 ) 



— 603234005661019847400a 14 c 3 - 41751662349615218000a 15 c 3 - 358769923659518800a 16 c 3 



968095828299600a 



17 + 



r (c_a+l) F ( c+a+50 ) 
,18 + 



32206016242400a i8 c 3 + 169856687000a iy c 3 + 106306200a 20 c 3 



+ - 



r (c_a+l) c+a+50 ) 

— 1401400a 21 c 3 - 2600a 22 c 3 + 138505441369671545827613807716230758400c 4 



F (c_a+1) r( C+a+50 ) 

— 141033120649691397441988498324822425600ac 4 



f c+a+50 



+- 



r(£=f±i) r( 2 , 

49191424454401897055805697928999731200a 2 c 4 

r(EfS) r( c+a + 50 ) 

— 7286784027262375583445490685229465600a 3 c 4 



- + 



F (C_C+ 1 ) F ( c+a+50 ) 



4„4 



374067588927556175588296387795353600a 4 c‘ 

r(£^|+i) r( c+a + 50 ) 



-+ 



+- 



12066879415392989959435569486336000a 5 c 4 - 1418367801582681846899974223616000a 6 c 4 

r( £ + ±i ) r( c+a 2 f50 ) 

— 11724885340790270723926442035200a T c 4 + 2225577070203302581076302156800a s c 4 

r(g=|+i) r( c+a + 50 ) + 

25892642522406895816692384000a 9 c 4 - 1660403587242510892624752000a 10 c 4 



r( 



c— a+1 



) r( 



c+a+50 



— 34718000105624284415366400a 11 c 4 + 349058892421463336112000a i2 c" 



12 „4 



r (c-a+l) c+a+50 ) 
,13 + 



- + 



15758070846864650304000a 13 c 4 + 98029228370495136000a 14 c 4 




28 



Salahuddin 



No. 1 



— 1569321907433280000a 15 c 4 - 25609443449990400a 16 c 4 - 91184052960000a 17 c 4 



523867344000a 18 c 4 + 4708704000a 19 c 4 



9609600a 2 °c 4 



47103137180828963726258689199082307584c 5 



-38606916894738607152732025212449587200ac s 



10872909845856838905945144160621363200a 2 c 5 



— 1250646162787845304931184049334476800a 3 c 5 

+ r(g=|+i) r( c+a + 50 ) + 

37484137237073610387082179113287680a 4 c 5 

+ r(£^|+i) r( c+a + 50 ) + 

2960056649806039350518792618176000a 5 c s - 147835927587805127993477304710400a 6 c 5 

r(^f^) r( c+a 2 f50 ) 

-4407777576011344380877953734400a T c 5 + 181587707284941430929852497280a 8 c 5 



5386052499736584815450220000a 9 c 5 - 73426510897677856161406800a 10 c 5 



— 3490081968585840976348000a 11 c 5 - 12978357558398102190320a 12 c 5 

+ r(£^|+i) r( c+a + 50 ) + 

7463477149491 16280000a 13 c 5 + 9937225737538642400a 14 c 5 + 2454621732763200a 15 c 5 

r(^f^) r( c + a + 50 ) 

— 738607265716320a ie c 5 - 4828187364000a 17 c 5 - 4859654800a ls c 5 + 39239200a 19 c 5 

f r(g=|+i) r( c + a + 50 ) 

i 80080a 2O c 5 + 11643512284338078214932291209134080000c 6 | 

+ r(£^|+i) r( c+a + 50 ) + 

— 7811846265289713988464106586033356800ac 6 



1782447335346527348520040105063219200a 2 c 6 



-156185991237696882904399849104998400a 3 c 6 



1628516051137267451572910992588800a 4 c b 



-8682292288286157632177578598400a 6 c 6 - 584337996155409991642835251200a 7 c 6 
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6839225908971453044632166400a 8 c 6 + 496205711149910054823129600a 9 c 6 



+ - 



r (c_a+l) r ( c+a+5U ) 

1041291387682609641177600a lo c 6 - 182911277098225554278400a n c 6 



r(c^a+l) r ( 



c+a+50 



-2197804342078164019200a i2 c 6 + 12956443561424563200a 13 c 6 



r(c = a+l) r( 



c+a+50 



-+ 



394556420664806400a 



2 / V 2 > 

14 - 6 1989973649664000a 15 c 6 - 6690587904000a 16 - 6 



r( 



c— a+1 



) r ( 



c+a+50 



-+ 



+- 



2 1 \ 2 > 

— 84756672000a 1 7c 6 - 192192000a 1 8c 6 + 2185434017754672893827554567467827200c 7 



r (c_a+l) r ( c+a+5U ) 

— 1210576372193996556349018510039449600ac 7 



r (c_a+l) r ( c+a+50 ) 



2233561 14720461583695372389167267840a 2 c 



2 „7 



F (c_a+ 1 ) F ( c+a+50 ) 



- + 



-14417840115170057077040994071756800a 3 c 7 - 127654808190213133895037229670400a 4 c 7 



c+a+50 ^ 



+- 



r(£-f±i)r( 2 

38362342001061342556663535539200a s c 7 — 100943627312800952414696545280a 6 c 7 



r (c_a+l) c+a+50 ) 

— 45930295596168058484716070400a 7 c 7 - 158519703458949712693632000a s c 7 



F (c = a+1) r( 



c+a+50 



+- 



26491220954206657796347200a 9 c 7 + 318516097490735070738240a lo c 7 



r( 



c— a+1 



)r( 



c+a+50 



-4796976148767599596800a n c 7 - 110735713090730041600a i2 c 



12 „7 



F (c_a+ 1 ) F ( c+a+50 ) 



+ - 



— 307118684131459200a 13 c 7 + 7303674051273600a 14 c 7 + 61439842464000a 15 c 7 



F (c_a+ 1 ) F ( c+a+50 ) 

86692320000a 16 c 7 - 504504000a 17 c 7 - 1144000a 18 c 7 



F (c_a+ 1 ) F ( c+a+50 ) 

320841683284681063430490974846976000c 8 - 147261685663349662978242993822105600ac s 



F (c_a+ 1 ) F ( c+a+50 ) 

21828588787720013764469751624499200a 2 c 8 - 978764458188092993647066703462400a 3 c 8 



F (c_a+1) F ( C+a+50 ) 

-29937782516869216684714136371200a 4 c 8 + 2477913960753858299794587648000a s c 8 



-+ 



r( 



c— a+1 



) r( 



c+a+50 



2 / \ 2 / 

30306732837796421151215616000a 6 c 8 - 2324406168836210117372313600a 7 c 8 



+ - 



F (c_a+ 1 ) F ( c+a+50 ) 

— 35079004867678815190118400a 8 c 8 + 810488727074359517184000a 9 c 8 



- + 



r( 



c— a+1 



) r ( 



c+a+50 
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18840273819769239552000a lo c 8 - 11075412804186931200a n c 8 

+ r(£^|+i) r (. c+a+50 ) + 

— 2870202835480780800a 12 c s - 20168019615744000a 13 c 8 + 36128581632000a 14 c 8 



774918144000a 15 c 8 + 1976832000a 16 c 8 + 37623854094699788409092847834234880c 9 



— 14310486246486860625946551595827200ac 9 + 1685406368191632966927765317222400a 2 c 9 



-46774053309260554179902883430400a 3 c 9 - 2907367898999056189584016384000a 4 c 9 

+ r(^=f^) r( c+a + 50 ) + 

110290157599761586095108096000a 5 c 9 + 3091438308324915618903142400a 6 c 9 

+ r(g=|+i) r( c+a + 50 ) + 

— 73000241928734665517772800a T c 9 - 2137975103648864229235200a 8 c 9 

+ r(^|+i) r( c+a + 50 ) + 

9101205323954551808000a 9 c 9 + 599787058295388620800a lo c 9 + 3364607652652646400a n c 9 

r(^f^) r( c+a + 50 ) 

— 36478880044006400a 12 c 9 - 421246938112000a 13 c 9 - 788646144000a u c 9 



3587584000a 15 c 9 + 9152000a 16 c 9 + 3578017643742734189270453452800000c 10 

r(^=f^) r( c + a + 50 ) 



-1124822208152756932629830094028800ac 10 + 103529588791248473877009294950400a 2 c 10 

r (e_a+l) c+a+50 ) 

— 1311215108521671070577904844800a 3 c 10 - 185621293383786411465336422400a 4 c 10 



3060077972270587621992038400a 5 c 1 



168536740537058850462105600a 6 c 1 



r (c_a+l) r ( c+a+50 ) 

— 1016343899030933407334400a 7 c 10 - 75771265876540543795200a 8 c 10 

r(c^|+l) c+a+50 ) 

— 312238666563644620800a 9 c 19 + 10967326380122112000a lo c 10 



110687000154931200a 11 c 1 ° - 58633012838400a 12 c 10 - 4098455961600a 13 c 1D 

r(c = |+i) r( c + a + 50 ) 

— 11948851200a 14 c 10 + 278986102350768311768576111411200c 11 



— 72124888712454702873593459507200ac 11 + 5065553589978581619224540610560a 2 c 1 



5329959949354438689384038400a 3 c 11 — 8535575626577387683612262400a 4 c 11 
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25100649148683283563315200a s c 11 



6059556534266084019322880a e c 11 



r (c_a+l) r( . e+a +5U ) 

22155491991978909491200a 7 c 11 - 1690530802659461734400a 8 c n 



r (£=f±i) r( c+a + 50 ) 



- + 



-16415314891981824000a 9 c n 



95043876381941760a 10 c n + 1706929285734400a 11 c 11 



- + 



r(c = a+l) r( 



c+a+50 ^ 



-,14 „11 



r( 



C— Q+l 



)T( 



c+a+50 



-+ 



2 7 V 2 > 

17973019590141218234189217792000c 12 - 3792934708034192465188198809600ac 12 

T+=fS) r( c+a + 50 ) 

196285879020447610271642419200a 2 c 12 



2 > " V 2 

2786214843661751921698406400a 3 c 12 



r (c_a+l) c+a+50 ) 

— 289801563760060600614912000a 4 c 12 — 2064705306471 178076160000a 5 c 12 



- + 



r(£=f±i) r( c+a + 5U ) 

148846508105105620992000a e c 12 + 1585317292830248140800a 7 c 12 



+ - 



r( £ +± i ) r( c+a + 50 ) 

— 22299454819830988800a®c 12 - 354397212868608000a 9 c 12 - 241573183488000a 10 c 12 



r (c_a+l) r ( c + a+50 ) 

45622886400a 12 c 12 + 961450433202482165846785392640c 13 



r (c_a+l) c+a+50 ) 

-163998227183988503762121523200ac 13 + 5924648410475117616444211200a 2 c 13 

r(^|+i) r( c+a 2 f50 ) 

166515545915840350896128000a 3 c 13 - 7213269573432725603942400a 4 c 13 
4 



+ - 



r (c_a+l) c+a+50 ) 

— 115037909882834792448000a 5 c 13 + 2404631605453956300800a e c 13 

r(c = |+i) r( c + a + 50 ) 



+- 



42663988630042214400a 7 c 13 - 110264137559654400a 8 c 13 - 4272497037312000a 9 c 13 



F (c_a+ 1 ) c+a+50 ) 



— 13576756838400a lo c 13 



- 42990796800a n c 13 



146227200a 12 c 13 



- + 



r (c_a+l) c+a+50 ) 

42819108327436116421509120000c 14 - 5828749054688624947023052800ac 14 



+ - 



r (c_a+l) F ( c+a+50 ) 

133998221660012686186905600a 2 c 14 + 6158133813035578687488000a 3 c 14 



- + 



r( 



c— a+1 



) r( 



c+a+50 



+- 



2 ) \ 2 > 

— 125789151018504683520000a 4 c 14 - 3283349863865057280000a 5 c 14 



F (c_a+ 1 ) F ( c+a+50 ) 



+ - 



21132740770489958400a 6 c 14 + 679691867258880000a 7 c 14 + 1409416298496000a 8 c 14 



r( 



c— a+1 



) r ( 



c+a+50 
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— 28005433344000a 9 c 14 - 114307891200a 10 c 14 + 1588225843907503003480883200c 15 



+ - 



r (c_a+l) r( . c+a + 5 U ) 

— 169762020930124137680076800ac 15 + 2043378481535587823124480a 2 c 15 



r (£=f±i) r( c+a + 50 ) 



162454135634000379904000a 3 c 15 - 1316763562211776921600a 4 c 



4„15 



r (c_a+l) r ( c+a+50 ) 

— 61099456493769523200a 5 c 15 - 21663400701788160a 6 c 15 + 6667626921984000aV 5 



r( 



C— Q+l 



)T( 



c+a+50 



2 / V 2 / 

28300728729600a 8 c 15 - 77792870400a 9 c 15 - 317521920a lo c 15 



+ - 



r (c_a+l) r ( c+a+ 5 U ) 

48967663640865914486784000c 16 - 4025992923028302436761600ac 16 

r(g=f+i) r( c+a + 50 ) 

12152046118550057779200a 2 c 16 + 3163650334517152972800a 3 c 16 



F (c_a+1) c+a+50 ) 



- + 



— 1616946504361574400a 4 c 16 - 767423554781184000a 5 c 16 - 3085932036096000a 6 c 16 



r (£=f±i) r( c+a + 50 ) 

37340577792000a r c 16 + 190513152000a 8 c 16 + 1249348776871366450216960c 17 



F (c_a+ 1 ) F ( c+a+50 ) 

— 76939510335577862963200ac 17 - 329033249485284966400a 2 c 17 



+ - 



F (c_a+1) F ( c+a+50 ) 

45434381571765043200a 3 c 17 + 199713870087782400a 4 c 17 - 6309093310464000a 5 c 17 



r( 



c— a+1 



) r ( 



c+a+50 



2 / v 2 / 

— 37687050240000a 6 c 17 + 91521024000a 7 c 17 



466944000a 8 . 



,17 



F (c_a+ 1 ) F ( c+a+50 ) 

26180456293639127040000c 18 - 1166186969396202700800ac 18 

+ F (g=|+i) r( c+a + 50 ) + 

— 11694085683491635200a 2 c 18 + 469496995774464000a 3 c 18 + 3535117221888000a 4 c 18 



F (c_a+ 1 ) F ( c+a+50 ) 

— 30751064064000a 5 c 18 - 209190912000a 6 c 18 + 445504624400741171200c 19 



r(£=f±i) r( c+a + 50 ) 

— 13687619181241958400ac 19 - 196960707910041600a 2 c 19 + 3308414369792000a 3 c 19 



F (c_a+1) F ( c+a+50 ) 

31005671424000a 4 c 19 - 67436544000a 5 c 19 - 458752000a 6 c 19 + 6054208717455360000c 20 



r( 



c— a+1 



) r( 



c+a+50 



2 / \ 2 / 

— 119879948147097600ac 20 - 2097585900748800a 2 c 2 ° + 14242598092800a 3 c 2 ° 



F (c_a+ 1 ) F ( c+a+50 ) 

145332633600a 4 c 20 + 64105204122910720c 21 - 737026192179200ac 21 



r( 



c— a+1 



) r ( 



c+a+50 
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-14349002342400 a 2 c 21 + 28259123200a 3 c 21 + 288358400a 4 c 21 + 509356277760000c 22 

r(sn|+i) r( c+a + 50 ) 

— 2836188364800ac 22 - 57881395200a 2 c 22 + 2855482163200c 23 - 5138022400ac 23 

+ r(g=|+i) r( c+a + 50 ) + 

— 104857600a 2 c 23 + 10066329600c 24 + 16777216c 25 
+ r(£=f±i) r(£±i|t50) 



( 8 ) 



§3. Derivation of summation formulae 

Derivation of main result (8): 

Substituting b = — a — 49, z = \ in given result (2), we get 



(2a + 49) 2-^1 



a , — a — 49 ; 1 

• 2 



= a 2 F\ 



a + 1 , — a — 49 ;1 

. 2 



49) 2 F\ 



a , — a — 48 ; 1 

. 2 



Now using same parallel method which is used in [6], we can prove the main result. 
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This class was introduced by Al- Amiri and Reade W. In particular Hq = R , the class of 
functions whose derivative has a positive real part and studied by Macgregor I 12 l . Also H\ = K. 

B a (a > 0) is the class of functions in A which satisfy the condition 

Re f'(z) (M) >0. (4) 

The class B a was introduced by Singh I 19 ! and studied further by Thomas I 21 ! and El- 
Ashwah and Thomas I 3 1. Functions of this class are called Bazilevic functions. Particularly 
B 0 = S* and B 1 = R. 

In 1976, Noonan and Thomas ^ 16 1 stated the qth Hankel determinant for q > 1 and n > 1 
as 

&n &n+l Q*n+q— 1 

«n+l 

O-n+q - 1 0‘n+2q-2 

This determinant has also been considered by several authors. For example, Noor 
determined the rate of growth of H q (n) as n — > oo for functions given by Eq.(l) with bounded 
boundary. Ehrenborg ^ studied the Hankel determinant of exponential polynomials and the 
Hankel transform of an integer sequence is defined and some of its properties discussed by 
Layman I 9 1 . Also Hankel determinant for different classes was studied by various authors 
including Hayrnan I 5 1, Pommerenke I 18 !, Janteng et alJ 6,7,8 ^ and recently by Mehrok and Singh 

[13] _ 

Easily, one can observe that the Fekete-Szego functional is ^(l). Fekete and Szego I 4 1 
then further generalised the estimate of |a 3 — /za| | where fi is real and / £ S. For our discussion 
in this paper, we consider the Hankel determinant in the case of q = 2 and n = 2, 

a-2 a 3 

«3 «4 

For fi complex, Szynal I 21 ! obtained the estimates for a , 3 — fj,a%\ for the class M a . Al- Amiri 
and Reade I 1 ! obtained the estimates for |a 3 — | for the class H a and Singh I 19 l obtained the 

estimates for |a 3 — fj,a% \ for the class B a . 

In this paper, we seek upper bound of the functional |a 2 a 4 — a|| for the functions belonging 
to the classes M ai H a and B a . Results due to various authors follow as special cases. 

§2. Main result 

Let P be the family of all functions p analytic in E for which Re(p(z)) > 0 [z £ E) and 

p{z) = 1 +PiZ +P 2 Z 2 + .... (5) 



H q (n ) = 



Lemma 2.1. If p € P, then \pk\ < 2 (k = 1, 2, 3, ...). 
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This result is due to Pommerenke I 18 ! . 
Lemma 2.2. If p £ P, then 



2p 2 =pI + (4 — Pi)x, 



4p 3 =p\ + 2pi(4 -p\)x -pi(k-p\)x 2 + 2(4 -p?)(l - \x\ ' 2 )z, 

for some x and z satisfying |x| < 1, \z\ < 1 and p\ € [0, 2]. 

This result was proved by Libera and Zlotkiewiez I 10 - 11 !. 

Theorem 2.1. If / £ M a , then 



|a 2 a 4 - ai\ < 



1 3a(l + a) 3 

(1 + 2a) 2 [ (1 + 3a) (2 + 15a + 24a 2 + 7a 3 ) 



+ 1 



Proof. As / € M a> so from (2) 



(l-«) 



zf\z) , _ (zf(z))' 



f(z) 



a- 



f\z) 



= p(z). 



On expanding and equating the coefficients of z, z 2 and z 3 in (7), we obtain 



a 2 



Pi 

1 + a 



= P2 (1 + 3 a)pl 

° 3 2(1 + 2a) ' 2(1 + 2a)(l + a) 2 

and 

P3 (1 + 5a)pip2 (4 + 6a + 17a 2 )p 3 

04 - 3(1 + 3a) + 2(l + a)(l + 2a)(l + 3a) + 6(1 + 2a)(l + 3a)(l + a) 3 ' 

Using (8), (9) and (10), it yields 



( 6 ) 



( 7 ) 

( 8 ) 

(9) 

( 10 ) 



|a 2 a 4 — c 



1 

C(a) 



4(1 + 2a) 2 (l + a) 3 p 4 (4p 3 ) + 12(1 + 2a)(l + 5a)(l + a) 2 p 2 (2p 2 ) 
+8(1 + 2a)(l + 6a + 17a 2 ) 2 pf — 3(1 + 3a)((l + a) 2 (2p 2 ) + 2(1 + 3 ce)p 2 ) 2 

(ID 



where C(a) = 48(1 + 3a)(l + 2a) 2 (l + a) 4 . 

Using lemma 2.1 and lemma 2.2 in (11) and replacing p 4 by p , it can be easily established 



that 



\a 2 a 4 -al\ < —— 
6(a) 



[-4(1 + 2a)(l + a) 2 - 12(1 + 2a)(l + 5a)(l + a) 2 - 8(1 + 2a)(l + 6a + 17a 2 )" 
+3(1 + 3a) (3 + 8a + a 2 )]/ + [8(1 + 2a) 2 (l + a) 3 + 12(1 + 2a) (1 + 5a) 
—6(1 + 3a)(3 + 8a + a 2 )(l + a) 2 ]p 2 (4 — p 2 )S + (1 + a) 3 (2 — p)[6(l + a)(l + 3a) 
— (1 + 4a + 7a 2 )](4 — p 2 )p 2 5 2 + 8(1 + 2a) 2 (l + a) 3 p(4 — p 2 ) 



where S = \x\ < 1. 

Therefore 

|a 2 a 4 - a|| < ^yP(<5), 



where 
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F(8) = [-4(1 + 2a) (1 + a) 2 - 12(1 + 2a)(l + 5a)(l + a) 2 - 8(1 + 2a)(l + 6a + 17a 2 ) 

+ 3(1 + 3a) (3 + 8a + a 2 )]p 4 + [8(1 + 2a) 2 (l + a) 3 + 12(1 + 2a)(l + 5a) - 6(1 + 3a) (3 + 8a 
+ a 2 )(l + a) 2 ]p 2 (4 — p 2 )8 + (1 + a) 3 (2 — ja) [6(1 + a)(l + 3a) — (1 + 4a + 7a 2 )](4 — p 2 )p 2 8 2 
+ 8(1 + 2a) 2 (l + a) 3 p(4 — p 2 ). 

As F'(6) > 0, so F(8) is an increasing function in [0, 1] and therefore MaxF(d) = F( 1). 
Consequently 

\a 2 a 4 -al\ < G(p), (12) 

C(a) 



where G(p) = F( 1). 

So 

G(P) 



— A(a)p 4 + B(a)p 2 + 48(1 + 3a)(l + a) 4 , 



where 



A(a) = 4a(l + a) (2 + 15a + 24a 2 + 7a 3 ) and B(a) = 48a(l + a) 4 . 

Now G (p) = —4 A(a)p 3 + 2 B(a)p and G (j>) = —12 A(a)p 2 + 2 B(a). 

G' (p) = 0 gives p[2A{a)p 2 — B(a )] = 0. 

G"(p) is negative at p = yj = P- 

So Ma xG(p) = G(p). 

Hence from (12), we obtain (6). 

The result is sharp for pi = p , p 2 = p\ — 2 and p^ = p\{2 — p\). 

For a = 0 and a = 1 respectively, theorem 2.1 gives the following results due to Janteng 
et alJ 8 J. 

Corollary 2.1. If f(z ) € S * , then 



| <220.4 — a 2 1 < 1. 



Corollary 2.2. If f(z) € K, then 



|a 2 a 4 - a 3 | < -. 



Theorem 2.2. If / G H a , then 



— a? I < < 9 ( 1+ “) 2 



|a 2 a 4 -o 3 |S 



(17a — 5) 2 



if 0 < a < yj. 
if jf < a < 1. 



144(l+2a)(l+20a+7a 2 -4a 3 ) ~ 9(1 +a) 2 

Proof. Since / € H a , so from (3) 

(1 - a)f(z) + = p(z). 

On expanding and equating the coefficients of z, z 2 and z 3 in (14), we obtain 

Pi 



a 2 = 



2 ’ 



P2 + api 
3(1 + a) 



(13) 



(14) 



(15) 



0-3 = 



(16) 
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and 



04 = 



P3 



3apip 2 



a(2a — 1 )p\ 



4(1 + 2a) 4(1 + a)(l + 2a) 4(1 + a)(l + 2a ) ' 

Using (15), (16) and (17), it yields 

9(1 + a) 3 pi(4p 3 ) + 54(1 + a)p\(2p 2 ) 
+36a(2 a - 1)(1 + a)pf - 8(1 + 2a) (2 p 2 + 2 ap\) 2 



|a 2 a 4 - a 3 | = 



1 

cW) 



(17) 



(18) 



where C(a) = 288(1 + 3a)(l + 2a)(l + a) 2 . 

Using lemma 2.1 and lemma 2.2 in (18) and replacing p\ by p, we obtain 



|a 2 a4-a 2 | < 



C(a) 



(1 - 12a + 3a 2 + 8a 3 )/ + 2(1 + 13a + 4a 2 )p 2 (4 - p 2 )S 
+(2 — _p) [16(1 + 2a) — (1 + 2a + 9a 2 )p 2 ](4 — p 2 )5 2 + 18(1 + a) 2 p( 4 — p 2 ) 



where 5 = \x\ < 1. 
Therefore 



where 



|a 2 a 4 - ai\ < 



C(a) 



m, 



F{8) = (1 - 12a + 3a 2 + 8a 3 )/ + 2(1 + 13a + 4a 2 )/(4 - p 2 )S + (2- 
p)[16(l + 2a) — (1 + 2a + 9a 2 )/](4 — p 2 )S 2 + 18(1 + a) 2 p{ 4 — /) 

is an increasing function. 

Therefore MaxF(d) = F{ 1). Consequently 



|a 2 a 4 - a 3 | < 



C(a) 



G(P), 



(19) 



where G(p) = F( 1). 

So 

G{p) = —A{a)p i + B{a)p 2 + 128(1 + 2a), 
where A(a) = 2(1 + 20a + 17a- 2 — 4a 3 ) and B(a ) = 4(1 + a) (17a — 5). 

Case I. For 0 < a < B(a) < 0. 

So G(p) is maximum at p = 0 and it follows the first result of (13). 

In this case, the result is sharp for pi = 0, p 2 = —2 and p 3 = 0. 

Case II. For / < a < 1, as in theorem 2.1, G(p) is maximum forp = y ( 4 +wa+ 7 a 2 - 4 a 3 j = 
p , so MaxG(p) = G{p ). 

In this case, the result is sharp for p 3 = p , p 2 = p\ — 2 and p 3 = p 3 (p 2 — 3) . 

Hence the theorem. For a = 0 in theorem 2.2, we obtain the following results due to 
Janteng et alJ 6 l. 

Corollary 2.3. If f(z) £ R , then 

o , 4 

|a 2 a 4 - a 3 | < -. 



Putting a = 1 in theorem 2.2, we get the following results due to Janteng et alJ 8 l. 
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Corollary 2.4. If f(z ) £ K, then 

\a 2 a 4 - oJ| < jjp 



On the same lines, we can easily prove the following theorem: 

Theorem 2.3. If / £ B a , then 



|o 2 a 4 - 03 1 < 



4 

(2 + a) 2 ' 



The result is sharp for p\ = 0, p 2 — —2 and P 3 = 0. 

For a = 0, theorem 2.3 gives the following result due to Janteng et alJ 8 l. 
Corollary 2.5. If f(z) £ S * , then 



1 0204 — 03! < 1. 

For a = 0, theorem 2.3 gives the following result due to Janteng et alJ 6 l 

Corollary 2.6. If f(z) £ R, then 

0 , 4 

|a 2 a 4 - o 3 | < -. 
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Abstract Mixed censored survival data arise from taking the maximum or minimum of two- 
component or even multi-component systems. Tied to this is the need to estimate maximal or 
minimal survival functions from related marginal survival functions. In this study we exploit 
the fact that the transformation of survival time dynamics of one component to another can be 
considered to result from variations in shift (mean), concentration (mode), spread (variance), 
degree of symmetry (skewness) and peakedness (kurtosis) as well as the underlying relative 
proportion of the censored data. Consequently, we propose some regression models that can 
serve for prediction purposes. The models are evaluated using a combination of some real data 
from the literature combined with some Weibull simulated data of varying levels of censoring 
and parameters. To the extent that the performance measures are excellent, the suggested 
models can be effective for prediction analysis. 

Keywords Maximal survival function, minimial survival function, mixed censored survival 
data, univariate. 
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§1. Introduction and preliminaries 

Mixed censored survival data arise from taking the maximum or minimum of two-component 
systems or even multi-components. Tied to this is the need to estimate maximal or minimal 
survival functions from related marginal survival functions. Bivariate survival data for two- 
component systems pervade many fields of scientific research. For example, medical researchers 
have obtained data on the infection times of the left and right kidneys of patients, and times 
to HIV infection of sex partners. Clayton (1978) applied regression model for censored survival 
data to epidemiological studies of chronic disease incidence. A related model for association in 
bivariate survivorship time distributions was proposed for the analysis of familial tendency in 
disease incidence. ^ The possible extension of the model to general multivariate survisorship 
distribution was indicated. Oakes (1982) discussed a reparameterization of a model introduced 
by Clayton in bivariate life-tables and the inference for the parameter governing the association 
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was considered when the marginal distributions are specified up to Lehnmann alternatives J 5 1 

The information matrix was derived explicitly and it was shown that the parameteriza- 
tion is moderately successful in introducing orthogonality between the associated parameters. 
Demographers have collected data on the two survival times of individual twins. Economists 
interested in modeling the joint retirement decisions of married couples have collected data of 
times to retirement for both the wife and the husband of each surveyed household. One of the 
common features shared by taking the maximum or minimum of these data is a mixture of 
the two duration times, calling for a joint distribution that allows a probabillity mass on the 
45° line in the sample space. The exploitation of facts of the transformation of survival time 
dynamics of one component to another which will result to variations in shift (mean), concen- 
tration (mode), spread (variance), degree of symmetry (skewness) and peakedness (kurtosis) 
as well as the underlying relative proportion of the censored data is the purpose of this study. 
Specifically, we propose some regression models that can serve for prediction purpose. 

In [4] , the likelihood function of the estimation from censored survival data was given and 
several examples are seen. The type of estimate studied here can be briefly indicated as follows. 
Suppose from a two-component system a random sample of N paired subjects are taken such 
that Th,T 12 , . . . ,Tl 7 v and T 2 i,T 22 , ■ ■ ■ ,T 2 n are two random variables. Taking the maximum 
and minimum between the observed lifetimes separately will result to two sets of mixed samples. 
The sample distribution functions Fy 2 ( t ) and Fy 1 (t) are naturally defined as that which assigns a 
probability of j? to each of the given values, so that Fy 2 ( t ) equals times the number of sample 
values less than the argument ty 2 and F Yl ( t ) equals ^ times the number of sample values less 
than the argument t Yl . Besides describing the sample, they are also non-parametric estimates 
of the mixed population distributions. When observations are incomplete the corresponding 
estimates are still step functions with discontinuities at the ages of observed deaths, but can 
no longer be obtained as a mere description of the sample. The two samples are incomplete in 
the sense that some items in the mixtures have given not random samples T 11; T 12 , . . . , Tjjv and 
T 2 i,T 22 , . . . ,T 2N of random variables T\ and T 2 itself (called the lifetimes), but the observed 
lifetimes Zi = (Zu,Z 2 i) where Zu = min(Tij, Cu) and Z 2 i = min(T 2 j, C 2 j); i = 1,2, 

Here Cu and C 2 i, called limits of observation (i.e. censored times), are constants or values of 
the other random variables, which are assumed to be independent of Tu and T 2i respectively 
(unless otherwise stated). For each item it is known whether one has 

T u < C u , Zu = Tu and T 2i < C 2 i, Z 2i = T 2i , (deaths), 
or 

Tu > Cu, Zu = Cu and X 2 i > C 2 i,Z 2 i = C 2 i, (losses). 

Ordinarily the Tu,T 2 i,Cu and C 2 i are so defined as to be necessarily non-negative. The 
mixed censored survival data is such that an item in the mixture belongs to one population with 
certain probability and to the other with another probability. This results into a division of 
mutually exclusive classes of deaths and losses. A loss by definition always precludes the desired 
knowledge of Ti = (Tu,T 2i ). On the other hand, a death does not preclude the knowledge of 
the corresponding C, = (Cu,C 2 i). 
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Notation 1.1. Let Xjj and Cu represent the failure and censoring time for the ith 
subject of the first component. Also let and Cu represent the failure and censoring time 
for the ith subject of the second component. Let the bivariate vector failure time T- t have 
cumulative distribution function Fj = ( Fu , X 2 j), probability density function /,; = (fu, / 2 i) and 
survivor function F t = (Xjj, Fu). Also let the bivariate censorsing life time Ci have cumulative 
distribution function, cummulative distribution function G* = (Gij,G 2 j), probability density 
function (ji = (gu,gu) and survivor function Gj = (Gij,G 2 j). Let Y be the corresponding 
vector of covariates. Let Tj = (Tjj,T 2 j) and C t = (Gij,G 2 j) be independently conditional on 
Yi . On each of N individuals, we observe the group (Yj, <5ij, 8 u, 1%) where = (Zu,Zu) = 
denotes the follow up time for the paired entity on individual i under study. 

Su = denotes a censoring indicator for the first member of an entity on individual i. 

5 2i = denotes a censoring indicator for the second member of an entity on individual i. 

7 i = the binary indicator used to indicate which of the following up time is considered on 
individual i. 

That is, 

f 1) if Z\ i < Z‘2i , 

H 

[ 0, if Zu > Z‘u- 

L m ax(Y, 61,62,1) = the likelihood function for the maximal survival function and L m i n (Y. , 61, 
62,”/) — the likelihood function for the minimal survival function. 

Observing the maximum between the follow up time for the paired entity on the ith individ- 
ual (i.e. Yi = max(Zu, Z 2 j), we have the likelihood function, for the N i.i.d groups (Y, 61,62,1), 
for the maximum survival function as 



<(Y, 61,62,1) = JJ{(5ii(y)-Fi*(y)) [(1 Al ‘ )(1 7,)] x (g2i(y)F 2 i(y)) [{1 <52,)7il 



Af2i(y)G2i(y)) [s ^ i] x {fiiWGuiv))^ 1 -^ 

[(1 — <52i)(l — 7*)] 



g2i{r)F 2i {T)dT 



/2i(T)G 2i (r)dr 



9u(T)Fii(T)dT 



l(l-«5ii)7i] 



[<52i(l-7i)l 



fii(r)Gu(T)dT 



[<5ii7i] 



}• ( 1 ) 



In the same light, observation of the minimum follow up time for the paired entity also on 
the ith individual (i.e. Yj = mm(Zu, Z21)) results in the likelihood function for the minimal 
survival function given by 



L m in(Y, 61,62, l) 



II {(9n(y)Fii(y)) [{1 - Slih,] x {92i{y)My)) l(1 - S2i)(1 - li)] 



i = 1 



x(f2i(y)Cu(y)) [s ^ 1 -^ x (fu(y)Gii(y))^ 

[(l-«2i)7il 



X / g 2 i(T)F 2 i{T)dT 



/ 2 i(r)G 2i (r)dr 



[<52.7 i. 



gii{T)Fu(T)dT 



fu(j )Gu(t )dr 



1(1-<5i.)(1-7.)] 



[<5ii(l-7i)[ 



}• ( 2 ) 
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Defined for various censored and failure situations. 

Notation 1.2. 



ho = 9i{y)Fi(y ) = Pr{Ci = y, T) > y}; * = 1,2, 

hi = fi(y)Gi{y) = Pr{Ti = y, Ci> y}\ * = 1,2. 
py _ roo 

Lij(y) = / hj(T)d,T and L^y) = / hj(T)dT\ i = 1, 2, j = 0, 1. 



J 0 J y 

It is clear that the terms of L max (Y, <$i, 5 2 , 7) and L m - ln (Y, Si, 62, 7) can be expressed as 



Max (Ml M 7) 



L mi n(Y, S 1 ,S 2 ,j) 



II{( / io(y)) [Cl - 5li)(1 - 7i)I x (i 2 0 (y)) K1 ~ 52ih ' ] 

i=l 

x(l 2 i(y)) [S2 ™ ] x (/ii(y)) [5li(1 - 7i)1 x (Li 0 (y)) [(1_<5li)7il 
x(L 2 o(j/)) [( 1 - 52 ' )(1 - 7i)I x (L 2 i(y)) [fe(1 - 7i)1 x (LnM)^]}, (3) 

n 

1 1 { (/ io(i/)) t(1 A;ihil x (/ 20 (y)) [(1 - 52i)(1 - 7i)1 

i=l 

X(l 2 i(y)f 2i(1 ^ i)] X (info))^ 1 "^ x (Z 10 (y)) [(1 “ 5li)7i] 
x(Z 20 (y))[( 1 -^^ 1 -^)] x x (in(!,))M}. (4) 



Dividing (3) by (4), we have (5) as 



L max (Y, 5 1 ,5 2 ,'y) 

Min(MiM7 ) 







[(1-i52)7] 



[<52(1-7)1 



L 



20 



(yA 



[(1 — <5 2 )(1— 7 )] 



^20 (2/) / 



-'ll 



(yA 



; n(y) / 

My) 

M(y) 



[<5i7] 



[<5 2 7] 



[10 (y) 
M(y) 
-^io(y) 
My) 



[(1-50(1-7)] 



[( 1 - 52 ) 7 ] ' 



(5) 



We introduce the following notations 



mij 



n 0 



My) 

hj(y) , 

-My)’ 



* = 1,2; j = 0,l, 

* = 1,2; 1=0,1. 



Remark 1.1. In light of this notations, likelihood functions here in this section, the 
ratio of the likelihood function of the the maximal and minimal survival functions will be 
given in Sections 2. It is established that any regression of maximal survival functions on 
minimal survival functions is essentially a relation of certain instantaneous censored and failure 
rates. The validation of these models (which will be seen in Section 2) was such that for 
fixed first component data of proportion p = of censoring to Weibull distribution (/?,?)), 
the second component data was generated for varied levels of censoring of proportion p = 
and || with mean parameters A and .75A. We thus generally have 6 combinations 
identified as (a, 6, c, d) = (mean of real first component data, proportion of censored points 
in first component data, mean of simulated second component data, proportion of censored 
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points in the second component data). The marginal, maximal and minimal survival function 
for the two components was obtained and appropriate regression analysis was carried out. The 
results will be presented also in Section 2 . Discussions and conclusions will be made. The 
study concludes that the eight models developed in this study gave good hypothesis and can 
reliably predict minimal survival function and maximal survival function from marginal survival 
functions and recommends that models such as the ones discussed herewith, in particular Models 
2 . 7 , 2.4 and 2.8 (the best performing models) be adopted for studying failure times in phenomena 
that can be modelled in terms of components. 

§2. Main results 

Using these notations in ( 5 ), we have the ratio of the likelihood of the maximal and minimal 
survival functions given by 



Ljnax{y, 8 l, 8 2 , 7 ) _ . 


f (m 10 (y)) (1 Sl ' > {rnii(y)) Sl \ 


1-7 
1 1 


f(m2o(y)) (1 52) {m 21 {y)) S2 \ 


L m in(Y,8i,8 2 ,7) 


V (^20 (y)) (1— (^21 (y) )^ 2 ) 


V («io(y)) (1_,5l) ( n ii(2/))' 51 J 



(6) breaks down, for specific values of 81,62,7 into the following ratio for likelihoods of 
maximum and minimal survival functions e.g. 

L max (y\8i =0, <$2 = 0, 7 = 0) _ rniojy) 

L m in(y\8i =0, 82 = 0, 7 = 0) n 20 (y) ' 

It is noted that rriij(y) and ntj(y), i = 1 , 2 ; j = 0, 1 can be likened to instantaneous 
event rates. (6) shows that the ratio of the likelihoods of the maximal and minimal survival 
function can be expressed as ratio of instantaneous censored and failure rates. Any regression 
of maximal survival functions on minimal survival functions is essentially a relation of certain 
instantaneous censored and failure rates. 

For the purpose of exploying this relation, it is necessary we state some results that suggest 
appropriate models of interest. 

Theorem 2 . 1 . Let Xi,X2 be stochastically independent random variables with distri- 
bution function Fi(-) and F 2 (-) respectively. Also, let maximum Y 2 = max.(Xi,X 2 ) and the 
minimum Y\ = min(Xi,X2) be ordered statistics respectively having the distribution func- 
tions G 2 (*) and Gi(-) respectively. Then for Y 2 and Y\ we have that G 2 (t) = Fi(t)F 2 {t) and 
G\(t) = 1 — ((1 — -Fi(£))(l — F 2 (t)). 

Proof. G 2 (t) = Pr(Y 2 < t) = Pr(X\ < t;X 2 < t) = Pr(Xi < t)Pr(X 2 < t ) (by 
stochastic independence) 

= F 1 (t)F 2 (t). ( 7 ) 

Similarly 

Gi(t) = 1 - Pr(Yi >t) = 1 - Pr{ X 1 > t, X 2 > t) = 1 - (1 - Fi(f))(l - F 2 {t)). (8) 

Model 2 . 1 . The (8) suggests the regression model In Fy 2i (t) = a± lnFu(t)+a 2 In F 2i (t) + 
In £i, where e* ~ !V( 0 , 1 ) and In = log e , F Yl (t) = Pr(Y 1 > t) = F m in( Xl ,x 2 ){t), Fx(t) = Pr(X 1 > 
t) and F 2 (t) = Pr{X 2 > t). 
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Model 2.2. The (7) suggests the regression model lnFy^i) = /3\ \nFu(t)+ f3 2 ln.F 2 j(f) + 
In Tj, where n ~ IV(0, 1) and In = log e , Fy 2 (f) = Pr(Y 2 > t) = F max{Xl ,x 2 )(t), F^t) = Pr(X 1 > 
t) and F 2 (t) = Pr(X 2 > t). 

Theorem 2.2. If F x {t) is a survivor function, then F Y {t) = > 0 is a 

proper survivor function for all 0 < /3 < 1. 

Proof. This follows trivially from the definition. 

From this relation, we have that 

1(1 ~F Y (t)) 1 ~Fx(t) 

F Y (t) F x (t ) ' 1 j 

(9) can be assumed to hold for 

1. F Y (t) = F max ( XliX2 ){t) against F x (t ) = F Xl (t), 

2. F Y (t) = F max{XuX2) {t) against F x (t) = F X2 (t), 

3. F Y (t) = F min(Xl ^ X2) (t) against F x (t) = F Xl (t ), 

4. F Y {t) = F min(XuX2) (t) against F x (t) = F X2 {t). 



In view of all these we formulate the following set of regression models. 

Model 2.3. 



1 

F Y 2i (t) 



7i- 



(1 ~ Fuji)) 
Fu(t) 



(1 ~F 2i (t)) . ^ 

72 F 2i (() +T ‘- 



where t* ~ iV(0, 1) and Fy 2 (f) = F max(Xl)X2) (t). 

Model 2.4. 



1 ~ ^Vii ft) 



Q!l 



(1 ~ Fuji)) 

F u {t) 



+ a 2 



(1 -_F Y2i (t)) 
F 2i {t) 



+ £*> 



where e* ~ A(0, 1) and Fy^f) = F min(Xl)X2) (t). 

Theorem 2.3. Assume that the random variable A is a symmetric triangular distribution 
with the probability density function symmetric about h given by 



h(x) = 



a-\-(x—h ) 



h — a < x < h, 



a—(x—h ) 



h < x < h + a, 



then the corresponding cumulative distribution is given by 

f 5^(2: -(/i- a)) 2 , 

H(x) = | 

1 2 ^(a 2 + 2a(x - h) - (x - /i) 2 ), 



h — a < x < h, 
h < x < h + a. 



Proof. See [1]. 

Theorem 2.4. Let fo(x) be a unimodal probability density function, FL{x) a cumulative 
distribution function having probability density function h(x), and assume both fo(x) and h(x) 
are symmetric about zero. 




Vol. 10 



Certain regression models involving univariate 
maximal and minimal survival functions on mixed censored survival data 



47 



Let f(x) = ^ fo( x )H(w(x)), — oo < x < oo and k = J x 2 fo(x)dx , then f(x) is a 

probability density function for any odd function w(x) and any a > 0. 

Proof. See [6]. 

Corollary 2.1. In the particular case that fo{x) in Theorem 2.4 is a probability density 
function w(x) = x, h( x) = |(1 — ^)i_ a a ( x ), so that 



The function 




2^2 (a 2 + 2 ax + x 2 ), 
2^2 (a 2 + 2 ax - x 2 ), 



—a < x < 0, 
0 < x < a. 



f{x) = 



al(i+afc) fo( x )( a2 + 2az + x 2 ), —a < £ < 0, 
a^a+afc) /oW(« 2 + 2ax - a: 2 ), 0 < x < a, 



is also a probability density function. I 1 ! 
Remark 2.1. We note as follows. 



1. the relationship in Corollary 2.1 can be simply expressed as f(x) = JT =1 a,iX l fo(x) for 
certain co-efficient a*, * = 1,2. 

2. The practical effect of Corollary 2.1 is to generate a possibly asymmetric and/or bimodal 
density distribution from a symmetric and unimodal density. 

3. A survival function may essentially be geometrically different from another survival func- 
tion by such characteristics (measures) as shift (mean), concentration (mode), spread 
(variance), asymmetry (skewness) and peakedness (kurtosis) all of which can be captured 
in the mechanism of the transformation in Models 2.5 and 2.6. 



Model 2.5. 

2 2 2 

ait l F u (t ) + ^2 bit l F 2 i(t ) + ^2 m it z Fu(t)F 2 i(t) + e i; where e* ~ N( 0, 1). 

i = 1 i = 1 i — 1 

Model 2.6. 

2 2 2 

fyiM = Cit l F u (t) + ^2 difF 2 i{t) + ^2 nit l F u (t)F 2i (t) + n, where n ~ N( 0, 1). 

2=1 i—1 i—1 

Theorem 2.5. If F(t) is a survival function, so also is, if F*(t) defined by 

/•-*(/) = r 1 ( ^<o )! ' 0 > 1. 

Proof. This follows from the definition of survival function. In view of Theorem 2.5 we 
postulate that P m ax(i) can be expressed in terms of F\ (t) and F 2 (t) as follows 
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Taking logarithm of (10), we have 

In F max (t) = 2 - (Fi(t)) 1 " 01 - {F 2 {t)Y~ e \ 

where 01,6*2 > 1 and In = log e . 

We can thus develop a regression model as such 

Model 2.7. 

In Fy 2i (t) = /3 0 +a 1 (F li (t))~i+a2(F li (t)y ; i+a3(F2i(t))~i+a i (F2i(t))~i+£i, where e* ~ 1V(0, 1). 

Model 2.8. 

In Fy u (t) = /3o+7i(-P , ii(6)) _ *+72(Ai(6)) _ *+73(^2i(i)) _ *+74(^2i(6)) _ ^+n, where n ~ N( 0,1). 

These sets of models were validated using a real data for first marginal component and 
simulated data for second marginal component. The parameters are generally somewhat cor- 
related, it may thus not be advisable to independently choose the parameters in a modelling 
situation. For example, the maximum likelihood estimate for the real data assuming Weibull ['1 
distribution for time to failure is given as /3 = 1.823823334 and fj = 5666.376110. In the light 
of these properties a set of data with some shape parameter (/?) but double the scale param- 
eter (77) will have double the mean. This point proved useful in the simulation exercise. The 
maximum likelihood estimates of the Weibull distribution for the real data were obtained and 
thus its mean (A, say). Considered to be of special interest is the extent to which these sets of 
models were validated to different mix (proportion) of censored data and for different closeness 
or balance of the marginal components. The proportion p of censored data were varied for 
different amounts p = 58/70,35/70 and 12/70. In the same way different combinations of the 
marginal components were considered. Thus, we had cases (A, A) and (A, 75A). This is to say 
that the distribution of the first marginal component is equal with that of the second marginal 
component in terms of the mean of the Weibull distribution and where there is a shift in the 
sense that the second marginal component is only .75 of the mean. This elucidated situations 
when censored points for the second component was high ( p = 58/70) just like the original 
data; censored points were even (p = 35/70) and censored points were sparse (p = 12/70). In 
summary, each of the regression models were validated with six generated samples of the sec- 
ond marginal component against the real data from the literature adopted for the first marginal 
component. Among the statistical tools that were employed as appropriate for the validation 
of the performances of these models were the adjusted R 2 (adj. Ft 2 ) and error sum of squares 
(SSE). The necessary results are given on Table 2.1 for various levels of censoring i.e. , highly, 
evenly and sparsely censored as well as for equal and shift in mean. Models 2.1, 2.3, 2.5 and 2.7 
connect maximal survival functions to the marginal survival functions while Models 2.2, 2.4, 

2.6 and 2.8 connect minimal survival functions to marginal survival functions. 

Table 2.1. Adjusted R 2 and Residual Sum of Squares Values for the Maximal and Mini- 
mal Response function by level of censoring of Simulated and Shift of Mean for the Models rela- 
tive to the real data (from combination (A, 58/70; A, 58/70) to combination (A, 58/70; .75A, 12/70). 
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Intensity (Mean) 


Response Variable 


Highly-censored 


Evenly-censored 


Sparsely-censored 






Adj. R 2 


SSE 


Adj. R 2 


SSE 


Adj. R 2 


SSE 


Equal Mean 


Model 2.1. 


.980 


.003 


.979 


.018 


.879 


.034 




Model 2.3. 


.977 


.004 


.964 


.054 


.968 


.191 




Model 2.5. 


.909 


.446 


.889 


.446 


.896 


.378 




Model 2.7. 


.992 


.000 


.997 


.001 


1.000 


.000 




Model 2.2. 


.985 


.003 


.951 


.018 


.998 


.010 




Model 2.4. 


.985 


.004 


.965 


.029 


.977 


.034 




Model 2.6. 


.904 


.470 


.901 


.468 


.912 


.391 




Model 2.8. 


.996 


.000 


.994 


.001 


.999 


.000 


Shifted Mean 


Model 2.1. 


.989 


.001 


.937 


.008 


.968 


.014 




Model 2.3. 


.992 


.001 


.979 


.011 


.971 


.019 




Model 2.5. 


.929 


.370 


.964 


.160 


.948 


.184 




Model 2.7. 


.993 


.000 


.963 


.004 


.935 


.009 




Model 2.2. 


.998 


.001 


.805 


27.814 


.993 


.083 




Model 2.4. 


.999 


.001 


.761 


1.844 


.999 


.333 




Model 2.6. 


.916 


.355 


.951 


1776.725 


.920 


.213 




Model 2.8. 


.999 


.000 


.998 


.003 


.989 


.071 



The study have been able to achieve the formulation of regression models based on the 
original thought that every survival function is a transformation of another survival function by 
such characteristics (measure) as shift (mean), concentration (mode), spread (variance), degree 
of symmetry (skewness) and peakedness (kurtosis). Based on the analysis of the data used, the 
following conclusions were reached. 

(1) all the sets of regression models can reliably predict minimal survival functions and max- 
imal survival functions from marginal survival functions. 

(2) irrespective of the level of censoring and shift in mean, Models 2.4. and 2.8. perform 
excellently in prediction of minimal response functions and the best overall is Model 2.8 
with adjusted R 2 in the range (.989 - .999). Model 2.8 compared with Model 2.6 in the 
range (.901 - .951). Similarly, in terms of residual sum of squares Model 2.8 is in the 
range (.000 - .071) and Model 2.6 is in the range (.213 - 1776.725) irrespective of whether 
it is equal or shifted mean. 

However for the prediction of maximal response function, the performance of Model 2.7 in 
relation to the others is less pronounced for shifted mean compared with when we have equal 
mean. This is especially so when using the criterion of the adjusted R 2 . It is noted that Models 
2.5 and 2.6 which makes allowance for bimodality have not performed well in relation to the 
others but it performed best where the proportion of censoring when the two components are 
closely equally censored (i.e. in the case where the data is approaching bimodal). 
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It is recommended that models such as the ones informed by (9) and (10), in particular 2.7, 
2.4 and 2.8 (the best performing models) be adopted for studying failure times in phenomena 
that can be modelled in terms of components. This is particularly so in the transformation 
when the survival time dynamics of one component to another component can be considered to 
be shift (mean), concentration (mode), spread (variance), degree of symmetry (skewness) and 
peakedness (kurtosis). 
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§1. Introduction 

By U, we denote the class of bounded analytic functions w(z) in the unit disc E = 
{z :| z |< 1} and of the form 

OO 

w(z) = ^2 c n z n , z G E, (1) 

n— 1 

which satisfying the conditions w{ 0) = 0, \w{z)\ < 1. 

It is well known that 

|ci|<l,|c 2 |<l-| Cl | 2 . (2) 

Let f(z) and F(z) be two analytic functions in the unit disc E, then f(z) is said to be 
subordinate to F(z) if there exists a function w(z) G U such that f(z) = F{w{z)) and we write 
as f{z) A F(z). 

In case the function F is univalent, the above subordination is equivalent to /( 0) = F(0) 
and f{E ) C F(E). 

Let A p (p is a positive integer) denote the class of functions of the form 

OO 

f(z) =z p + J2 an +P z n+p , ( 3 ) 

n— 1 

which are analytic in the unit disc E. 

By A, we denote the class of functions of the form 

OO 

f{z) = z + '22a n z n , 

n—2 



( 4 ) 
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analytic in the unit disc E. 



S* represents the class of p-valently starlike functions f(z) in A p and satisfying the condi- 



tion 



Re 



zf(z) 

[ f(z) J 



> 0, z G E. 



S*(A,B) represents the class of functions f{z) in A. p which satisfy the condition 

zf{z) -< pAAl, -i<b<a<i,zge. 



/W 



1 + Bz’ 



( 5 ) 



( 6 ) 



Clearly S*(l, —1) = S*, which shows that S*(A, B) is a subclass of S*. To avoid repetition, 
we lay down, once for all that — 1 < B < A < 1 and z £ E. 

S* is the class of analytic univalent starlike functions f(z) in A such that 



Re 



zf(z) 

m 



> 0. 



S*(A,B) is the subclass of univalent starlike functions f(z ) in A such that 

zf'(z) t 1 + Az 
f(z) * 1 + Bz‘ 



( 7 ) 



( 8 ) 



Janowski I 6 1 introduced the class S*(A,B) and proved some results by the method of 
variational techniques. Goel and Mehrok I 2,3 ! studied the class S* ( A , B) using subordination 
principle. 

The class K p consists of functions f(z) in A p for which 



Re 



(zf'(z)Y 

f'(z) J 



> 0. 



(9) 



Clearly f(z) € K p implies that z ^ p z ^ £ S*. 

K p (A , B) is the subclass of functions f(z) in A p for which 

(zf(z)y l + Az 
f'(z) P l + Bz 

K stands for the class of convex univalent functions f(z) in A with the condition 

'( zf(z)r 



Re 



f'(z) 



> 0. 



K (A, B) is the subclass of K of functions f(z) in A with the condition 

(zf(z)y ^ 1 +Az 



f'(z) 



1 + Bz 



( 10 ) 



(ID 



( 12 ) 



Let M p {a){a > 0) denote the class of functions f(z) in A p which satisfy the conditions 



± 0 and 



Re 



zf(z) (. zf(z)y 

(1 — a) — — b a 



f(z) 



f'(z) J 



> 0. 



( 13 ) 
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M p (a; A, B)(a > 0) denotes the subclass of M p (a) of functions f(z) in A p and satisfying 



the conditions Y 1 0 and 



„ ^zf(z) , Jzf(z)Y , J + Az 

(1 — a ) \ + a JTT\ < P-i , n ’ 
f{z) f'(z) 1 + Bz 



(14) 



M(a ) denotes the class of functions f(z) in A and satisfying the conditions ^ Y 0 and 

(15) 



Re 



M ,zf(z) , (*/'(*))' 

{1 ~ a) lW + a ^W 



> 0. 



The class M (a) was introduced by Mocanu I 9 1 . Miller et al. have shown that all a-convex 
functions are starlike in E and for a > 1, all a-convex functions are convex in E. 

M(a;A,B)(a > 0) is the subclass of M(a) of functions f(z ) in A and satisfying the 
conditions ^ Y 0 and 



(1 — a ) 



zf(z) , (zf(z)Y , J + Az 
f{z) f'(z) <P l + Bz- 



(16) 



The class M(oe;A,B) was studied by Goel and Mehrok I 4,5 !. Fekete and Szego I 1 ! made an 
early study for the bounds of |a 3 — /xa| | (/x real) when f(z) is analytic univalent in E. The well 
known result due to them states that if f(z) = z + a nZ n is analytic and univalent in E, 

then 



( 3 — 4/i. i_i < 0; 



\a 3 - na\\ < < 



1 + 2 exp{j^), <)<//< 1 : 



[ 4/i — 3, fi > 1. 



The following observations are obvious: 

(i) M p (a; 1 , — 1 ) = M p (a). 

(ii) M 1 (a-,A,B) = M{a-,A,B). 

(hi) M p (0-A,B) = S;(A,B). 

(iv) M p (0; 1 , — 1 ) = S*. 

(v) M 1 (0-A,B)=S*(A,B). 

(vi) Mi(0; 1,-1) = S*. 

(vii) M p (l-A,B) = K p (A,B). 

(viii) M p ( 1 ; 1 , - 1 ) = K p . 

(ix) M 1 (1-A,B) = K(A,B). 

(x) Mi(l ;!,-!) = K. 
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( 17 ) 

( 18 ) 

( 19 ) 

( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 

( 24 ) 
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Using (2) in (24), we obtain 



\a p+2 ~ /XGp+ll < 



(A - B)p 2 (A - B) 2 p 4 

2(p + 2a) + ( p + a ) 2 



[\\- p\ -v]. 



where v is defined in (19). 

Case I. /i < A. 

From (2.9), we have 



2 (A — B)p 2 (A — B) 2 p 4 

K +2 - pa p+1 < 2{p+2a) + {p + a)2 [(A - *) - A*]. 



If /x < (A — v), from (26) 



- 



If A — v < p < A, again from (26) 



i 2 i ^ ( A - B )P 2 

|a p+2 pa p+ 1\< 2{p + 2a y 



Case II. p > A. 

From (2.9), we have 



2 . (A — B)p 2 ( A — B) 2 p 4 

|a p+2 - /m p+1 | < 2(p + 2a) + (p + a)2 \M - (A + *)]• 



If /x < (A + z/), from (27) 



, 2 ^ (A- B)p 2 

\a p+2 pa p+1 \ < 2(p + 2Q) • 



If M > (A + z/), again from (27) 

2 , ( A-B) 2 p 4 . 

|«p+ 2 - Ma p+ i| < ( p + a ) 2 (A* - A )- 

For a = 0, equality sign hold in first and third inequalities of (17) for the function fi(z) 
defined by 



fi(z) = 



z p (l +B6 z) E ^ M1 , B^ 0; 



{ z P e (P ASz \ B = 0 , \S\ = 1 . 

and equality sign hold in second inequality of (17) for the function fo{z) defined by 

f z p {l + BSz 2 ) si ^ir 1 , By£ 0 ; 
fo{z) = < 

z p e s *^ 1 ,B = Q,\6\ = 1. 
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For a > 0, equality sign hold in first and third inequalities of (17) for the function / 1 (a)( 2 ) 
defined by 

t)] a ,B^ 0 ; 



/l(a) i z ) 



s(./o tS Hi + B5t) Pl B~o. B> dtj 
[h(fo^- 1 e iS ^ 1 ctt)] a ,B = 0,\S\ = l. 

and equality sign hold in second inequality of (17) for the function /o( Q )(~) defined by 



(30) 



/o(a)(H 



(fjt* 1 (1 + BSt 2 ) P< 2Bc > dt 



,B? 0; 



1 e pA ^ ) dtj , B = 0, |<5| = 1. 



(31) 



Remark 2.1. Choosing p = 1 in the Theorem 2.1, we get the results proved by Goel and 
Mehrok M. 

Choosing A = 1,B = — lin the Theorem 2.1, we have 
Corollary 2.1. If / £ M p (a), then 



\ a p+2~ MOp+ll < { 



where 



and 



Ai — 



(^p(A 1 — At), M < Ai - zq; 
2 

1 Al — Vl < p < Ai + V\\ 

(p+a) 2 - Al), p > Ai + I/!. 

2p(p 2 + 2po + a) + (p + a) 2 
4p 2 (p + 2a) 

(p + a) 2 

Vl = 



(32) 



4p 2 (p + 2a) ' 

For a = 0, the extremal function for the first and third inequalities of (32) is obtained 
by putting A = 1 } B = ~ lin (28) and extremal function for the second inequality of (30) is 
obtained by putting A = 1, B = — 1 in (29). 

For a > 0, the extremal function for the first and third inequalities of (32) is obtained 
by putting A = 1, B = — lin (30) and extremal function for the second inequality of (32) is 
obtained by putting A = 1, B = — 1 in (31). 

Choosing a = 0 in the Theorem 2.1, we have 

Corollary 2.2. If f £ S*(A,B), then 



( p(A-B) 
2 



U - B)p -B]-p 2 (A~ B) 2 p , p < {A - 2 f A p _-^ Bh , 



\a p+2 - pa 2 p+1 \ < < 



p(A-B) (A-B)p-( 1+B) < < (A-B)p+(1-B) . 

2 ’ 2 (A-B)p — B — 2 (A-B)p ’ 



(33) 



[ P 2 (A ~ B) 2 p - — B)p — B], p > (A ~ 2 3- + i)r~ a) • 

The first and third inequalities of (33) are sharp for the function fi(z) defined in (28) and 
the second inequality is sharp for the function fo(z) defined in (29). 
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Remark 2.2. By taking p = 1 in the Corollary 2.2, we get the estimates of |a 3 — pa\ \ 
for the class S*(A, B). 

Remark 2.3. Taking A = 1, B = —1 in the Corollary 2.2, we obtain the following result: 
If f £ S*, then 

( p(2p+l) - 4p 2 /q n<\] 



\a p+2 - juap +1 | < 



P, \ < P- < 



p+i . 

2 p ' 



{ \jrp-p(‘2p- n> . 



Remark 2.4. Taking p = 1, A = 1, B = — 1 in the Corollary 2.2, we obtain the following 
result: If f € S*, then 



3 - 4/z, ix < §; 



|«3 - pal\ < < 



l, b <n< l; 



[ 4/i — 3, fi > 1. 



This result is due to Keogh and Merkes I 7 ! . 
Putting a = 1 in the Theorem 2.1, we get 
Corollary 2.3. If / £ K p (A,B), then 



( (. A-B)[(A-B)p-B]p 2 _ (A-B)W u < (r+l) 2 [(A-B)p-(l+B)] , 

j 2(p+2) (p+1) 2 ^ - 2(A-B)p 2 (p+2) ’ 



l a p+ 2 M a p+ll — 



(A-B)p 2 (p+l) 2 [(A-B)p-(B+l)] < , (p+l) 2 [(A-B)p+(l-B)] . 

2(p+2) > 2(A-B)p 2 (p+2) - t 1 — 2(A— B)p 2 (p+2) ’ 



(34) 



(,4-B)V (A-B)[(A-B)p-B]p 2 ^ (p+1) 2 [(A— B)p+(1 — B)] 

k (p+1) 2 B 2(p+2) ’ B — 2(A— B)p 2 (p+2) 



The first and third inequalities of (34) are sharp for the function defined by 



h{z) = 



f 0 z ptP-^l + B8t)‘^B m dt, B/ 0; 



f z pt p ~ 1 e pASt dt, B = 0, |<5| = 1. 



and the second inequality of (34) is sharp for the function defined by 



( f~ pt p 7 (1 + B5t 2 ) F '' is ) dt, B/0; 

fo(z) = l 

[ f~pt p ~ 1 e pA * t dt , B = 0, |<5| = 1. 



Remark 2.5. Letting p = 1 in Corollary 2.3, we obtain the estimates of |«3 — pa 2 \ for 
the class K(A, B). 

Remark 2.6. 




58 



B. S. Mehrok and Gagandeep Singh 



No. 1 



Letting A = 1, B = — 1 in Corollary 2.3, we have the result: If f £ K p , then 



\a p+2 - pa 2 p+l \ < < 



p 2 (2p+l) 4 p 4 ^ (p+1) 2 . 

p+2 (p+l) 2 ^’ ^ — 2p(p+2)’ 

P 2 (P+1) 2 < .. < (P+1) 3 ■ 

p+2 7 2p(p+2) - - 2p 2 (p+2)> 



4p 4 „ P 2 (2p+1) > (p+1) 3 

(p+1) 2 1 1 p+2 > " — 2p 2 (p+2) ’ 



Remark 2.7. Taking p = 1, A = 1, B = — 1 in the Corollary 2.3, we obtain the following 
result: If f € K, then 

1 - H, M < |; 



l«3 -M^l < | < M< t 

. M - 1 , M > f 



3’ 



This result is due to Keogh and Merkes . 
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over cube-full number 
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Abstract Let n > 1 be an integer, q^\n) is a multiplicative function. In this paper, we shall 
study the mean value of exponential divisor function qf\n) over cube-full number, that is 

( lk\ n ) = ^W/sCn), 

n<x n<x 

n is cube — full 

where f 3 (n) is the characteristic function of cube-full integers, i.e. 

{ 1, n is cube-full; 

0, otherwise. 

Keywords cube-full number, divisor problem, Dirichlet convolution, residue theorem, mean 
value. 



§1. Introduction 

Let n > 1 be an integer of canonical form n = p^p^ 2 ■ ■ • p“ r . The integer n is called a fc-full 
number if n = pfp^ 2 ■ ■ -p“ r , where a\ > k,a ,2 > k, ■ ■ ■ ,a r > k. Let fk(n) be the characteristic 
function of fc-full integers, i.e. 

{ 1, n is fc-full; 

0, otherwise. 



We call an integer n = nl_i exponentially fc-free if all the exponents aj(l < i < r) are 
fc-free, i.e. are not divisible by the fc-th power of any prime (fc > 2). Let qjf* (n) denote the 
characteristic function of exponentially fc-free integers. 

Obviously the function q^\n) is multiplicative, and for every prime power p a , there are 
<lk\p ) = Qk\p 2 ) = 9k\p 3 ) = ■■■ = ?i e) (P 2fc_1 ) = l ^k\p 2k ) = 



^^This work is supported by Natural Science Foundation of China (Grant No: 11001 154), and Natural Science 
Foundation of Shandong Province(Nos:BS2009SF018, ZR2010AQ009). 
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Many authors have studied the properties of the exponential divisor function (ri). L. 
Toth proved the following result: 

^2 q[ e \n) = D k x + 0{x 1/2 5(x)), 

n<.x 

where 

P a= 2 k y 

qk(n) denoting the characteristic function of fc-free integers. 

In this paper, we shall study the mean value of exponential divisor function fjf (n) over 
cube-full number, that is 

E «i e) (n) = 

n<.x n<.x 

n is cube— full 

where fs (n) is the characteristic function of cube- full integers, i.e. 



h(ji) 



1, n is cube-full; 
0, otherwise. 



Theorem 1.1. For some D > 0, 



E 

n<x 

n is cube— full 



M, n) c(|)C(!)G(|) * , c(|)C(!)G(|) i , c(|)C(|)G(i) i 

% \ n ) — rvor X 3 



c(!) 



C(2) 



- x 4 



c(f) 



X 5 



+ O (x s exp ( — D(log x) t (log log x) ® ) ^ , 



where G(s) = ft (l — ^57 H ) is absolutely convergent for Res > 

p 

Notations 1.1. Throughout this paper, e always denotes a fixed but sufficiently small 
positive constant. 



§2. Proof of the theorem 

In order to prove our theorem, we need the following lemmas. 

Lemma 2.1. i 3 l 

A(3, 4, 5; x) <C log 3 x. 

Proof. The proof of this bound depends on the theory of two-dimensional exponent 
pairs. 

Lemma 2.2. Suppose f(n) is arithmetical function, and satisfy 

i 

Y J f{n) = Y. xaip A logx) + 0(x“), 

n<x j = 1 

E i/( n )i = °( xai lo s ra 0> 

n<.x 
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over cube-full number 



here ai > <12 > ••• > a; > 1/c > a > 0, r > 0, Pi (£),••• , p(t) are polynomials on t whose 
degree are not exceed r, and c > 1, b > 1 are fixed integers. 

Suppose for Res > 1, have 



y- Hb(n) 

2-^i n s 

n = 1 



1 

C WY 



if 



h(n) = Y Vb{d)f(n/d c ), 

d c | n 



then 



1 

Y' J /i(n) = ^ x aj Rj (log x) + E c (x), 

n<x j = 1 

here Ri(t), ■ ■ ■ , Ri(t) are polynomials on t whose degree are not exceed r, and for some D > 0, 
have 



E c {x) < £ 1/c exp (( - L»(loga:) 3/5 (logloga;) _1/5 )). 



Proof. If b = 1, see theorem 14.2 in A. Ivic I 2 1. If b > 2, it can be proved in the same 
method. 

Lemma 2.3. Let s = a + it is a complex number, then we have 



E 

n—1 

n is cube— full 



dk\n) 

n s 



C(3s)C(4s)C(5s) 

C(8s) 



G(s), 



( 1 ) 



where the Diriclrlet series G(s) = II (l L + ■ ■ ■ ) is absolutely convergent for Res > |. 

P 

Proof. Let 

00 

F(s) := E 

n—1 
n is cube 



n. s 2 -^ n.s 



— full 



n—1 



here 



h{n) 



1, n is cube-full; 
0, otherwise. 
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From q^\n) is multiplicative we have for Res > 1 that, 
F(s) = 



E 

n— 1 
cube- 

nfi+ 



Qk\n) = -y q { k\ri)h(n ) 

rr 1 S ' rv - > S 



n— 1 

n is cube— full 



n — 1 






+ 



IT p3s p4s p5s 



pis pis p. 



p2s 


+ - 


b 


1 ^ 


p(2 fc — l)s y 


1 


1 


p2 k s 


p(2 fc +l)s 



n3s 



+ ' '• + 



p(2 fc — l)s 



C(3s)C(4s) n(l + ^ 5s p 2 fc s p( 2 fc H 



(2 fc +l)s 



= C(3s)C(4s)C(5s lIK'-i-i 

p F ^ 

_ C(3s)C(4s)C(5s) tt /.. 1_ \ 

C(8a) 111 



C(3a)C(4 S )C(5 S ) 

C(8s) 



G(s), 



where G(s) =11(1“ H ), an< 4 it is absolutely convergent for a > 4 + e. 

p 

Next we prove our theorem. 

From lemma 2.3, we have known 

F(s) :=£?£! M = C(3»)a4»X(5,) G(8)i 



n— 1 

n is cube— full 



C(8s) 



where G(s) is absolutely convergent for a > | + e. 
Define 



G( S ) := £ 



s(«) 



n— 1 



oo , / \ oo X] d(3,4,5;m)5(0 

i?(s) := C(3s)C(4s)C(5s)G(s) := ^ Er = , ( Res > 1) 

z — ' n s z — ' n s 

n = 1 n=l 

here /i(ra) = X] d(3, 4, 5; m)g(l). Then we can get 



i=ml 



E M») = E ^(3,4,5 ;m)g{l) = ^2g(l) ^ rf(3,4,5;m). 



n<tc ml<x _ 

From Residue theorem and lemma 2.1, we can get 



E rf (3,4,5;n) = C(|)C(|)^ +C(|)C(j)®^ + <(|)C(^M + A(3,4,5;x) 



n<x 



^.4.^.5. i s/3, 5. i .,3, .,4, i 22 . 

= C( ^ 3 )' E 3 +C(|)C(4)® 4 +0(* m ). 



( 2 ) 
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Then from (2) and Abel integral formula, we have the relation: 

5>(n) = £s( 0 [c(^)C(f)(f)* + C(f)C(^)(f)* + C(f)C(^)(f)* +o((f)^) 

n<a? Kx 



= «!)«!)** £ $ + C( |)C( >» £ ® f + 4® £ It^) 



5(0 



,3, 5 , 



5(0 



3, _,4, 



5(0 



15(01 



Kx 

oo 



Kx 



Kx 



Kx 



= «!)«!)*» £ + C(|)«|M £ f§ + C(|)C(|)-* £ f§ + ■ o(xS £ Ifffl) 

/=1 Z=1 Z=1 / >tc 

+ °0' £ w) + °0' £ £ 

Z>x Z># Z<a? 



15(01 

[22/177 



= C(^K(^)G(i)xt + C(!)C(|)G(|)a:* + C(|)C(|)Gf(|)** + 0(xi). 

From lemma 2.2 and Perron’s formula M, we can get: 

) f „\ _ CCMM1),I - c(S)C(S)g(j) ,» ^ c(§)c(f)gq) ,i 



5fc( n ) = 



n<x 

n is cube— full 



C(l) ~ ' C(2) C(f) 

+ o(x* exp ( — D(loga;)s(logloga:) _ 5)^ . 

Here D > 0, and G(s) = JJ (l ^ ^ H ) is absolutely convergent for Res > |. 

p 

Now our theorem is proved. 



x 5 
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Abstract The Smarandache function S(n) is defined as the minimal positive integer m such 
that n|m!. The main purpose of this paper is to study the analyze converges questions for 

oo oo 

some series of the form , i.e., we proved the series s(n) s diverges for any 5 < 1, 

n=l n = 1 

oo 

and converges for any e > 0. 

n= 1 ' ' 

Keywords Smarandache function, smarandache series, converges. 



§1. Introduction and results 

For every positive integer n, let S(n) be the minimal positive integer m such that n|m!, 

i.e., 

5(n) = min {m : in £ N, n\m\}. 

This function is known as Smarandache function I 1 !. Easily, one has 5(1) = 1,5(2) = 2,5(3) = 
3, 5(4) = 4, 5(5) = 5, 5(6) = 3, 5(7) = 7, 5(8) = 4, 5(9) = 6, 5(10) = 5, • • • . 

Use the standard factorization of n = p^p^ 2 • • • p^ k , pi < P 2 < • • • < Pk , it’s trivial to 

have 

5(n) = m|K {5(p“‘)}. 

Many scholars have studied the properties of 5(n), for example, M. Farris and P. Mitchell 
I 2 1 show the boundary of S(p a ) as 

(p - l)a + 1 < S(p a ) < (p— l)[a + 1 + log p a] + 1. 

Z. Xu I 3 1 noticed the following interesting relationship formula 




by the fact that S(p) = p for p prime and 5 (n) < n except for the case n = 4 and n = p, where 
n(x) denotes the number of prime up to x, and [x] the greatest integer less or equal to x. Those 
and many other interesting results on Smarandache function 5(n), readers may refer to [2]- [6]. 
1 This work is supported by by the Science Foundation of Shanxi province (2013JM1016). 
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Let p be a fixed prime and n € N, the primitive numbers of power p, denoted by S p (n), is 
defined by 

S p (n) = min {m : m € N,p”|ro!} = S{p n ). 

OO 

Z. Xu obtained the identity between Riemann zeta function ((s) = -^,<7 > 1 and an 

n—1 

inhnite series involving S p {n) as 

y- 1 _ CO) 

Spin) P s ~ 1’ 

and he also obtained some other asymptotic formulae for S p (n). F. Luca I 4 1 proved the series 

OO OO 

* 5 converges for all S > 1 and diverges for all <5 < 1, and the series ^ s( 

n= l S '' n > U n - 1 

converges for any £ > 0. 

In this note, we studied the analyze converges problems for the inhnite series involving 
S(n). That is, we shall prove the following conclusions: 

Theorem 1.1. For any S < 1, the series 



E 



i 



-1 S(nY 



diverges. 

Theorem 1.2. For any e > 0, the series 



E 



S(n) eS ^ 



converges. 



§2. Some lemmas 



To complete the proof of theorems, we need two Lemmas. 

Lemma 2.1. Let p be any fixed prime. Then for any real number x > 1, we have the 
asymptotic formula: 



E 

n—l 

S p (n)<x 



l 

Sp{n) 



1 

P~ 1 



^ln x + 7 + 



plnp\ 

p-l) 



+ 0{x~i + e), 



where 7 is the Euler constant, e denotes any fixed positive numbers. 

Proof. See Theorem 2 of [3]. 

Lemma 2.2. 1 7 1 Let e > 0 and d{n) denotes the divisor function of positive integer n. 
Then 

d(n) = 0{n e ) < C e n e , 
where the o-constant C e depends on e. 

Proof. The proof follows [7] by writing n = Y\p a , the standard factorization of n. Then 

p\n 



p ae > 2 ae = e aeln2 > ae In 2 > ^(a+ l)eln2. 




66 



Qianli Yang 



No. 1 



If p e > 2, then p ae > 2 a > a + 1. Therefore, 



d(n) -J— r CX + 1 "I r cx + 1 -j— r CX ~b 1 ^ -p-r 

XX pOte XX pae XX pae — XX 

p\n p\n p\n p\n 



V i(a+l)eln2 V a + 1 

p\n z p\n 



p e < 2 p e >2 



The last item in above inequality is which is less than pp-^ = C e , say, the 



o-constant C e depends on e. 



e In 2 1 10 11 e In 2 

p|n p e <2 

p e < 2 



§3. Proof of theorems 



Proof of Theorem 1. 



We may treat the case 5 = 1 first. By Lemma 1 and the notation S p (n) = S(p n ), we have 



OO 1 OO 1 

E l li \ ^ I 

Q ( 7l n \ t. — ^--1-00 • J 



^ S(p n ) 



tc^+oo S v (n) 

n — 1 y ' 



S p (n)<x 



Obviously, for 6 < 1, Y s(n) d ' di ver g es follows easily by the trivial inequality: 

n= 1 

OO ^ OO ^ OO ^ 

^ S(n) s - ^ S{n j - ^ S(p n ) ’ 

n = 1 ' n—1 v 7 n—1 

complete the proof. 

Proof of Theorem 2. 

OO 

It certainly suffices to assume that e < 1. We rewrite series Y as 



where u(k) = |j {n : S(n) = k}. For every positive integer n such that S(n) = k is a divisor of 
k\, i.e. u(k) < d(k\). By Lemma 2 and the inequality bellow 






we have 



u{k) < d{k\) < C e (k\) e < C e 



where C e means that the constant depending on e. 

Therefore, recalling that the properties of the sequence (l + j :) k , we have 

^u(k) i fk + lY k 1 fk + l\ ek „ ^ 1 



C_XJ / 7 \ (-XJ 






( -' 1 X/ 9£fe ’ 
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for some constant Ci, it follows that series is bounded above by 

k=i 



completing the proof. 




fc= 1 
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Abstract A generalization of the Weibull geometric distribution so-called the transmuted 
Weibull geometric distribution is proposed and studied. Various structural properties includ- 
ing explicit expressions for the moments, moment generating function of the new distribution 
are derived. The estimation of the model parameters is performed by maximum likelihood 
method. We hope that the new distribution proposed here will serve as an alternative model 
to the other models which are available in the literature for modeling positive real data in 
many areas. 
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§1. Introduction and motivation 

The Weibull distribution is one of the most commonly used lifetime distribution in model- 
ing lifetime data. In practice, it has been shown to be very flexible in modeling various types of 
lifetime distributions with monotone failure rates but it is not useful for modeling the bathtub 
shaped and the unimodal failure rates which are common in reliability and biological studies. 
Several distributions have been proposed in the literature to model lifetime data. Adamidis and 
Loukas introduced the two-parameter exponential-geometric ( EG ) distribution with decreas- 
ing failure rate. Kus ^ introduced the exponential-Poisson distribution (following the same idea 
of the EG distribution) with decreasing failure rate and discussed various of its properties. Mar- 
shall and Olkin ^ ' 1 presented a method for adding a parameter to a family of distributions with 
application to the exponential and Weibull families. Adamidis et al. I 1 ! proposed the extended 
exponential-geometric ( EEG ) distribution which generalizes the EG distribution and discussed 
various of its statistical properties along with its reliability features. The hazard function of 
the EEG distribution can be monotone decreasing, increasing or constant. Barreto-Souza et 
al. introduced a Weibull-geometric (WG) distribution following from compounded a Weibull 
distribution with a geometric distribution and obtained a new lifetime distribution with more 
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flexibility than the Weibull distribution and study some of its properties. Suppose that {Yi}f =1 
are independent and identically distributed ( iid ) random variables following the Weibull dis- 
tribution W(a,9) with scale parameter a > 0, shape parameter 9 > 0 and probability density 
function 

g(x\a,9) = 0ot e x 8 - 1 e~( ax ' > ,x > 0, 

and N a discrete random variable having a geometric distribution with probability function 
P(n;p) = (1 — p)p n ~ 1 for n G N and p G (0, 1). Let X = min The marginal pdf of X 

is given by 

f{x;p, a, 9) = 9a 6 (1 — p)x e ~ 1 e~ ( - ax ' ) 1— pe~^ ax ^ ,x>0, (1) 

which defines the WG distribution (see Barreto-Souza et al. I 6 1. It is evident that (1) is 
much more flexible than the Weibull distribution. The EG distribution is a special case of 
the WG distribution for 9=1. When p approaches zero, the WG distribution leads to the 
Weibull W(a,6) distribution. The corresponding cumulative distribution function and hazard 
rate function are given by 

1 — e~^ ax ^ 

F(x-,p,a,9)= ——^, x >0, (2) 

and 

h(x-,p,a,9) = 9a e x e ~ 1 1 — , (3) 

respectively. 

In this article we defined the family of transmuted Weibull-geometric distribution. The 
main feature of this model is that a transmuted parameter is introduced in the subject distri- 
bution which provides greater flexibility in the form of new distributions. Using the quadratic 
rank transmutation map studied by Shaw et al. , we develop the four parameter transmuted 
Weibull-geometric TWG(x;p,a,9, A). We provide a comprehensive description of mathemati- 
cal properties of the subject distribution with the hope that it will attract wider applications 
in reliability, engineering and in other areas of research. The concept of transmuted explained 
in the following subsection. 

§2. Transmutation map 

In this subsection we demonstrate transmuted probability distribution. Let F\ and F 2 be 
the cumulative distribution functions, of two distributions with a common sample space. The 
general rank transmutation as given in Shaw et al. I 18 l is defined as 

Gri 2 {u) = F 2 {F{ 1 (u)) and G R2 i(u) = Fi(F 2 _1 (m)). 

Note that the inverse cumulative distribution function also known as quantile function is 
defined as 



F 1 (y)= infxeR {^(a;) > y} for y G [0, 1] . 
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The functions Gr 12 (w) and Gr 2 \(u) both map the unit interval I = [0,1] into itself, and 



under suitable assumptions are mutual inverses and they satisfy Gmj( 0) 
A quadratic Rank Transmutation Map (QRTM) is defined as 


= 0 and G Rij ( 0) = 1. 


Gri 2 (u) = u + Xu(1 - u), A < 1, 


(4) 


from which it follows that the cdf’s satisfy the relationship 




F 2 (x) = (1 + X)F 1 (x) — XFi(x) 2 , 


(5) 


which on differentiation yields, 




f 2 (x) = f 1 (x) {(l + X)-2XF 1 (x)} 


(6) 



where fi(x) and f 2 (x) are the corresponding pdfs associated with cdf F\(x) and F 2 {x) respec- 
tively. An extensive information about the quadratic rank transmutation map is given in Shaw 
et al. I 18 l that at A = 0 we have the distribution of the base random variable. 

Many authors dealing with the generalization of some well- known distributions. Aryal and 
Tsokos M defined the transmuted generalized extreme value distribution and they studied some 
basic mathematical characteristics of transmuted Gumbel probability distribution and it has 
been observed that the transmuted Gumbel can be used to model climate data. Also Aryal and 
Tsokos PI presented a new generalization of Weibull distribution called the transmuted Weibull 
distribution . Recently, Aryal I 18 l proposed and studied the various structural properties of the 
transmuted Log- Logistic distribution, Muhammad Khan and King P°1 introduced the trans- 
muted modified Weibull distribution which extends recent development on transmuted Weibull 
distribution by Aryal et al. PI and they studied the mathematical properties and maximum 
likelihood estimation of the unknown parameters. Merovci introduced the transmuted Rayleigh 
distribution, transmuted generalized Rayleigh distribution, and transmuted Lindley distribu- 
tion and they studied the mathematical properties and maximum likelihood estimation of the 
unknown parameters. Elbatal PI introduced transmuted modified inverse Weibull distribution. 



§3. Transmuted Weibull geometric distribution 



In this section we studied the transmuted Weibull-geometric (TWG) distribution. Now 
using (1) and (2) we have the cdf of transmuted Weibull-geometric (TWG) distribution 



^ TWG (*G Pi ^ 1 



1 _ e -(oix) 9 

1 -pe-( ax '> e 



1 + A - A 



1 - 

1 — pe~( ax ) e 



(7) 



where A is the transmuted parameter. The corresponding probability density function (pdf) of 
the transmuted Weibull-geometric is given by 



A) = /wg(x) [(1 + A) - 2A F wg (x)} 



= to 6 (l-p)i 8 - 1 e»‘ 



1 — pe 



-(o 



(1 + A) — 2A 



t \ & 

1 - e~^ 

1 — pe~( ax ^ e 



x 



( 8 ) 
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respectively. 

Figures (1) and (2) illustrates some of the possible shapes of the pdf and cdf of the TWG 
distribution for selected values of the parameters 9,p, a and A, respectively. 





Figure 1: The pdf’s of various TWG distributions for values of parameters:fig a) 9 = 
l;a = 1.5A = 0.5, fig b) 9 = 2.5, 3, 3.5, 4, 4.5, 5, 5.5; A = 0.5, 0.6, 0.7, 0.8, 0.9, 0.4, 0.3 and for 
p = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, a = 1.5; with color shapes purple, blue, plum, pink, red, black and 
darkcyan, respectively. 





Figure 2: The cdf’s of various TWG distributions for values of parameters: :fig a) 9 = 

l;,a = 1.5A = 0.5, fig b) 9 = 2.5, 3, 3.5, 4, 4.5, 5, 5.5; A = 0.5, 0.6, 0.7, 0.8, 0.9, 0.4, 0.3 and for 
p = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, a = 1.5; with color shapes purple, blue, plum, pink, red, black and 
darkcyan, respectively. 
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The reliability function (RF) of the transmuted Weibull-geometric distribution is denoted 
by Rtwg{x ) also known as the survivor function and is defined as 



Rtwg{x ) = 1 - F T wg{x) 

1 - 

l-pe-( aa d 9 



1 + A - A 



1 - \ 

1 — pe~( ax ') 6 J 



( 9 ) 



It is important to note that Rtwg{x) + Ftwg{x) = 1. One of the characteristic in 
reliability analysis is the hazard rate function (HF) defined by 



h-TWG(x) 



frwGjx) 

1 - Ftwg(x) 



( 10 ) 



It is important to note that the units for h,TWG{x) is the probability of failure per unit of 
time, distance or cycles. These failure rates are defined with different choices of parameters. 
The cumulative hazard function of the transmuted Weibull-geometric distribution is denoted 
by Htwg(x) and is defined as 



Htwg{x) = - In 



/ \ 0 

1 _ e -(«x) 

1 -pe-( ax ') e 



1 + A - A 



/ i _ e -(«*)» \ 

y 1 — pe~^ ax '> e J 



( 11 ) 



It is important to note that the units for Htwg(x) is the cumulative probability of failure 
per unit of time, distance or cycles, we can show that. For all choice of parameters the 
distribution has the decreasing patterns of cumulative instantaneous failure rates. 

Figure 3 illustrates some of the possible shapes of the hazard function of the TWG distri- 
bution for selected values of the parameters 6,p , a and A, respectively. 





Figure 3: The hazard functions of various TWG distributions for values of parameters:fig a 
9 = 1; , a = 1.5A = 0.5, fig b: 9 = 2.5, 3, 3.5, 4, 4.5, 5, 5.5; A = 0.5, 0.6, 0.7, 0.8, 0.9, 0.4, 0.3 and 
for p = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, , a = 1.5; with color shapes purple, blue, plum, pink, red, black 
and darkcyan, respectively. 
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§4. Statistical properties 



This section is devoted to studying statistical properties of the ( TWG ) distribution. 

4.1 Moments 

In this subsection we discuss the r t h. moment for (TWG) distribution. Moments are nec- 
essary and important in any statistical analysis, especially in applications. It can be used to 
study the most important features and characteristics of a distribution (e.g., tendency, disper- 
sion, skewness and kurtosis). 

Theorem 3.1. If X has TWG(x\ </>),</>= ( p , a, 9 , A) then the r t h moment of X is given 
by the following 



MrO) = (! + 

3=0 



(j + i)(i- P ym + i) 



"(. i + 



- A Z 

3=0 



(i-p)m + i) 



U + 2) (3 + 1) 



(J + 1 ) 9 U + 2 ) t 



( 12 ) 



Proof. Let X be a random variable with density function (8) . The rth ordinary moment 
of the (TWG) distribution is given by 



H r (x) = E(X r > = / x 1 f(x,<f>)dx 



= 9a e (l~p) < (1 + A) / x r+9 ~ 1 e~ < - ax ' >e 



I —2 



oo 

- 2A J ^+0-i e -(«*) 9 (l- e -(“ x)e ) (l -pe- (ax)9S ) 



-3 



If \z\ < 1 and k > 0, we have the series representation 

(1 r(A; + j) j 

(1 ^ Z) 

Expanding the term ^1 — pe~^ ax ^ ^ as in (14) yields 



{ OO 00 

(l + A )^(j + iy J x r+o-i e -U+iKax) 9 dx 

j=o 0 

00 °? 

- A J2(j + 1)0' + 2 V / x r+9 - 1 e- {j+1)(ax '> 9 (1 - e- (ax)9 ) 

a — n J 



~ (j + i)(i- P yr(! + i) ^ 

- (1 + A) z^ ^ (j + i)? 



dx 

J 

(1 -P) r (§ + 1) 



i=o 



l=o 



a' 



' O' + 2) _ (J + 1) ' _ 
.U + !) S (j + 2)£_ ’ 



(13) 



(14) 



( 15 ) 




74 



F. Merovci and I. Elbatal 



No. 1 



Which completes the proof. 

Based on the first four moments of the (TWG) distribution, the measures of skewness 
7L($) and kurtosis /c($) of the (TWG) distribution can obtained as 

_ /+3(0)-3/ii(0)/i 2 (g) + 2/i?(0) 

and 

= /i 4 (g) - 4/zi(fl)/z 3 (6>) + 6 nl(0)n 2 (0) - 3/4(6*) 

^ [/i 2 (g)-/i?(g)f 

4.2 Moment Generating function 

In this subsection we derived the moment generating function of (TWG) distribution. 
Theorem 3.2. If X has (TWG) distribution, then the moment generating function 
Mx(t) has the following form 



r— 0 



Mx(t) = E)d V 1 + A ^E 

3=0 



(i + i)(i- p ym + i) 



^ A E 

3=0 



(i- P )m + i) 



a r (j + l)e 

(j + 2) (j + 1) 



(16) 



. (j + 1) ® (j + 2) 0 J 

Proof. We start with the well known definition of the moment generating function given 



by 



00 

M x (t) = J e tx fTWG(x,4>)dx 



OO -Lf OO .ip 

= E X xr fTWG(x, (f>)dx = /4 



r — 0 



r — 0 



= E^a + A)E 

3 = 0 



r— 0 

(j + i)(i- P ym + i) 



- A E 

3=0 



(i-p)r(g + 1 ) 



ye 

a r (j + 1)5 

(j + 2) (j + 1) 



(j + 1) 9 (j + 2) ^ J 



(17) 



Which completes the proof. 



§5. Distribution of the order statistics 

In this section, we derive closed form expressions for the pdfs of the rth order statistic of the 
(TWG) distribution, also, the measures of skewness and kurtosis of the distribution of the r t h 
order statistic in a sample of size n for different choices of n; r are presented in this section. Let 
Xl,X 2 , . . . ,X n be a simple random sample from (TWG) distribution with pdf and cdf given 
by (7) and (8), respectively. Let Xl,X 2 , . . . ,X n denote the order statistics obtained from this 
sample. We now give the probability density function of X r:n , say f r - n (x, $) and the moments 




